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ABSTRACT 

We present photometry and derived redshifts from up to eleven bandpasses for 9927 galaxies in 
the Hubble Ultra Deep field (UDF), covering an observed wavelength range from the near-ultraviolet 
(NUV) to the near-infrared (NIR) with Hubble Space Telescope observations. Our Wide Field Camera 
3 (WFC3)/UV F225W, F275W, and F336W image mosaics from the ultra-violet UDF (UVUDF) 
imaging campaign are newly calibrated to correct for charge transfer inefficiency, and use new dark 
calibrations to minimize background gradients and pattern noise. Our NIR WFC3/IR image mosaics 
combine the imaging from the UDF09 and UDF12 campaigns with CANDELS data to provide NIR 
coverage for the entire UDF field of view. We use aperture-matched point-spread function corrected 
photometry to measure photometric redshifts in the UDF, sampling both the Lyman break and Balmer 
break of galaxies at z ~ 0.8 — 3.4, and one of the breaks over the rest of the redshift range. Our 
comparison of these results with a compilation of robust spectroscopic redshifts shows an improvement 
in the galaxy photometric redshifts by a factor of two in scatter and a factor three in outlier fraction 
over previous UDF catalogs. The inclusion of the new NUV data is responsible for a factor of two 
decrease in the outlier fraction compared to redshifts determined from only the optical and NIR data, 
and improves the scatter at z < 0.5 and at z > 2. The panchromatic coverage of the UDF from the 
NUV through the NIR yields robust photometric redshifts of the UDF, with the lowest outlier fraction 
available. 

Subject headings: cosmology: observations — galaxies: high-redshift — galaxies: photometry - 
galaxies: distances and redshifts — galaxies: evolution — 
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1. INTRODUCTION 


The Hubble Ultra Deep Field (UDF; iBeckwith et al.1 
is one of the most studied fields on the sky, 
with extremely sensitive high-resolution imaging cov¬ 
ering many photometric bandpasses. While only 
covering 12 arcmin 2 in the sky, the data’s depth, 
resolution, and wavelength coverage enable a wide 
range of scientific work. The data have been used to 
measure the colors a nd luminosity funct i on of high 
redsh i ft galaxies (e.g.. IBouwens et all 12006 1 IRvan et all 


20071 ICameron fe Drived 20091 Bouwens et al.l 2011 


S 


Dunlop et al.1 120 lT Schenker et all 20131 McLure et al l 

2013HQesch et al.ll2013l : IBouwens et al.ll20 1 -11 1 , to observe 


the evolution in the star form ation rate density (e.g., 
IBouwens et al.l l2007. 2009 j2010lh to dete rmine the mor- 


2004 lElmegreen et al. 

20054 Straushn et al.l 

2006; 

Hathi et al. 2008a]; 

Elmegreen & Elmegreen 

2010; 

Oesch et al.l 2010bt lElmeexeen & Elmesreenl 20141), to 


constrain the s t ar formation efficien c y of gas at z ~ 3 
dWolfe fc Chenl 120061 : iRafelski et al.l feoilll . to charac- 
terize new types of galaxies, such as c lumpy galaxies 
(|Elmeereen et al.1 120071 : IGuo et akl 12012ft . and for many 
other investigations. 

The upgrade of the Hubble Space Telescope (HST) 
in the fourth servicing mission added the Wide Field 
Camera 3 (WFC3) to HST’s instrumentation, en¬ 
abling very deep near-infrared (NIR) (UDF09 and 
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UDF12; lOesch et al.l l2010dlbl: IBouwens et al.l l2011bl : 
iKoekemoer et al.1 1201 .'1 lEllis et al. 20 131) and near- 
ultraviolet (NUV) fUVUDF: ITeplitz et al.l I2013T) ob¬ 
servations of the UDF. The NIR data enable the 

study of galaxies at the h ighe st redshifts at z _ > 

7 (e.g., lOesch et al.1 [2010dlbl: IBouwens et, al l 

Finkelstein et al.1 120101: IBouwens et al.1 l2011bl 


2010 ; 


2013 


Ellis et al.l 120131 : iDunlop et al.l I2013H . while the NUV 


observations enable studies of the rest-frame ult ravio¬ 
let (UV) at intermediate redshifts, z ~ 1 dTe plitz_et_al 


20Kj: lBond et al.ll20l4lKurczvnski et al.ll20l4lMei et al 


2014). 


Reliable redshifts of galaxies in the UDF are needed for 
such studies, but only a small number of galaxies in the 
UDF have spectroscopic redshifts, beca use of the faint¬ 
ness of the galaxies sampled (see Section ITTl) . Therefore, 
most studies in the UDF rely on color selection tech¬ 
niques or photometric redshifts, or use a small number 
of galaxies. While color selection techniques are very 
useful for sele cting a specific typ e of galaxy in a redshift 
interval (e.g.. iSteidel et al]l2003lh photometric redshifts 
have the advantage of determining redshifts for a large 
sample of galaxies, making use of all the photometric 
bandpasses simultaneously, while prov iding uncertain¬ 
ties on the resultant redshifts (e.g.. iBenftezlfeOOCil ). The 
only public photometric redsh ift catal og c overing the en¬ 
tire UDF is presented in I Coe et all (|2006f l, and was cre¬ 
ated before the new NIR and NUV observations of the 
UDF. Since then, other redshift catalogs for subsets of 
the data have been released, but cover fractions of the 
field of view (FOV) and smaller numbers o f galaxies (e.g., 
iRafelski et al .1120091 : iCameron et al.l[20Tlh . 

Photometric redshifts are best determined when in¬ 
cluding strong features in the spectral energy distribu¬ 
tions (SEDs), and galaxies exhibit multiple such features 
enabling robust redshift determinations. The most ap¬ 
parent feature is the Lyman break, composed of the Ly¬ 
man limit at 9I2A and the Lyman series lines short- 
ward of 1216A. This is followed by the Balmer break at 
3646A, and then the 4000A break, which is composed 
of the Call H and K lines and the s udden onset o f pho- 
tospheric opacity by ionized metals (lHamiltonlll985li . 

This paper presents photometric redshifts that include 
high resolution NUV and NIR data, complementing the 
optical data and thereby covering the wavelength range 
0-2 < A < 1.8/zm, providing the best dataset to measure 
redshifts in the UDF. The data enable measurements of 
both the Lyman break and the Balmer breaks simulta¬ 
neously in the interval 0.8 < z < 3.4. Additionally, 
at least one of these breaks, along with a long base¬ 
line measurement of the SED slope, is observed over 
the entire redshift range. This effectively removes the 
redshift degeneracies observed when limited to optical 
data, and thereby significantly reduces the outlier frac¬ 
tion when comparing photo metric redshifts t o a spectro¬ 
scopic redshift sample (e.p ., iCoe et aLll200f& lllbert et,~ahl 
2006; Rafelski et a 11 [2001 


These are the first photometric redshifts that make use 
of the new WFC3 NUV data, which improve the redshifts 
by sampling the Lyman break of high redshift galaxies, 
and more clearly defines the 4000A break for low red¬ 
shift galaxies. The NUV photometry requires careful 
calibration to obtain robust photometry, due to chal¬ 


lenges with charge transfer efficiency and dark calibra¬ 
tions. The methods to overcome these challenges and 
produce the high quality NUV image mosaics, which are 
now available on MAST0 are described, and the NUV 
photometry is provided. 

In addition to the NUV data, we include NIR data cov¬ 
ering the entire UDF. Previous work included the deep 
NIR data from the UDF09 and UDF12 campaigns; in 
this paper the photometric redshifts make use of both 
of these deep NIR data and shallower NIR observations 
covering the entire UDF FOV from CANDELS as de¬ 
scribed below. The NIR data measures the 4000A break 
at intermediate redshifts and provides more photometric 
points to constrain the SED slope, thereby improving the 
photometric redshifts. 

The inclusion of these two new datasets make the red¬ 
shifts estimated here the most accurate and robust red¬ 
shifts available at this time, with significant improve¬ 
ments for most of the ten thousand galaxies in the UDF. 
The photometry and redshifts tabul ated here have al¬ 
ready been used by various studies dBond et, al.1 12014 
iKurczvnski et al.l 12014 iMei et all 120141 . and are now 
available to the public. 

The paper is organized as follows: In Section [2l we de¬ 
scribe the observations, we explain the new calibrations 
of the NUV data, and we characterize the data. In Sec¬ 
tion [3] we detail the methodology for aperture-matched 
PSF photometry and catalog definition. In Section d] 
we determine the photometric redshifts, and assess their 
quality and improvements with the addition of the NUV 
data. We describe the final photometric and redshift cat¬ 
alog in Section [5j and summarize the paper in Section [Gl 
The Appendix includes a careful description of the need 
for specialized dark calibrations for WFC3/UVIS and the 
methodology used to implement them. 

2. OBSERVATIONS 

This study makes use of eleven photometric band- 
passes covering the UDF (a(J2000) = 03 ,t 32 m 39 s , 
d(J2000) = —27°47'29."1) at high spatial resolution, 
spanning wavelengths from the NUV to the NIR. The 
blue box in Figure [l] outlines the NUV coverage of the 
UDF (UVUDF), which is comprised of three WFC3- 
UV IS filters: F225W, F 275W. and F336W, as described 
by (ITeplitz et al.l 1 20 lit ). These data h a ve be en modified 
from those discussed in ITeplitz et al.l (1201311 to include 
improved calibrations, photometry, and astrometry as 
described below. The optical data defines the field region 
in the Figure, and is covered by the four original ACS 
optical filters: F435W . F606W, F775W, and F850LP 
(iBeckwith et all 12006 ). The red box in the Figure out¬ 
lines the deep NIR coverage, which includes four WFC3- 
IR filters: F105W, F125W, F140W, and F160 W ob¬ 
tained in the UDF09 and UDF12 program s (Ocsc h~et al.1 
2010dlbl: IBouwens et, al.l 1201114 IKoekemoer et, al.l 120131 : 
Ellis et, all 1201311 . The entire field is also covered by 
three of the four WFC3-IR filters (F105W, F125W, 
and F160W) in the CANDELS GOODS-S o bservations 
(iGroein et al.1 120111 : IKoekemoer et al.l I2011H . The full 
UDF field, NUV coverage, and deep NIR coverage cover 
areas of 12.8 arcmin 2 , 7.3 arcmin 2 , and 4.6 arcnrin 2 re¬ 
spectively. The full UDF coverage area shrinks to 11.4 

18 http://archive.stsci.edu/prepds/uvudf/ 
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Figure 1. UDF coverage maps for the WFC3-UVIS, ACS- 
WFC, and WFC3-IR data overlaid on a color image created 
from all eleven bandpasses. The entire field is image d in 
the original four ACS optical filters JBeckwi th et al .] [200ql and 
the three WFC3-IR fil ters from CANDELS dGrogmeHahl 120111 ; 
IKoekem oer et all 1201 11). The UV coverage for the three WFC- 
UVIS filters lTeplitz et all 1120131 1 is outlined by the blue rectangles, 
and the deep IR regi on for the four WF C 3-IR filters is outlined 
by the red rectangle lO esch et alil201 0dlbt IBouwens et al.M2011bl ; 
IKoekemoer et aLl 2013: Ellis et al. 20lSb The green rectangle out¬ 
lines the area covered by more than 30% exposure time in the 
optical UDF. North points upwards in the image. 


arcmin 2 when considering only the area covered by more 
than 30% of the optical exposure time, as shown by the 
green box in the Figure. 

We do not use the NIR NICMOS J and H data 
([Thompson et al.l l2005ll or the WFPC2 F300W obser- 
vations dVover et al.1 1200911 in our analysis because the 
same wavelengths are covered by other filters presented 
here with superior depth, resolution, and coverage area. 
Also, in general we do not include lower resolution 
observati ons of the UDF, suc h as ground-based imag¬ 
ing fe.g.. [Rafclski et al.l 120091 ) or low-resol ution space- 
based imaging such as Spitzer IRAC (e.g.. lAshbv et al.l 
l2013lb to avoid source confusion, correlated photometry, 
and the added systematic uncertainties associated with 
more complex photometric te chniques (see Sect ion 13.11) . 
Lastly, the SBC far-UV data dSiana et al.ll2007T ) are also 
not used here as they contain very few sources and thus 
would not significantly affect the results. The filter cov¬ 
erage and throughput for all eleven bandpasses is shown 
in Figure [2] The large number of filters over such a large 
wavelength range enables the simultaneous measurement 
of the Lyman Break and the Balmer/4000A break at red- 
shifts 0.8 < z < 3.4, significantly improving the photo¬ 
metric redshifts of these galaxies (see section H. 

2.1. General UDF Imaging Information 

The data represent the deepest high-resolution 
panchromatic dataset available, although the properties 
of each bandpass varies by instrument and filter. Table 
|T] provides information about each bandpass, including 
the effective wavelength, zero point, Galactic extinction, 
exposure time, depth, areal coverage, and point spread 



Wavelength (A) 


Figure 2. Total system throughput of the WFC3-UVIS, ACS- 
WFC, and WFC3-IR filters used in the UDF. Each filter is plotted 
as a separate color, with lines of the same color denoting the wave¬ 
length range sampled. The gray bars at the top show the wave¬ 
length range used for each of the three instruments. Note that 
these throughputs include the quantum efficiency of the CCD. 

function (PSF) full width at half maximum (FWHM). 
The Galactic extinction is de r ived from the extinction 
coefficients in iPostman et al.l (120121). using EfB—Vi = 
0.00782 based on the iSchlegel et all (1998) IR dust- 
emission maps. The PSF FWHM is measured by fit¬ 
ting a sy mmet rical Gaussian to the PSFs as described 
in Section 12.71 Although the PSFs are not Gaussian, 
these provide estimates of the resolution of the images 
for reference purposes. 

The depth measurements in Tabic [T| are obtained by 
measuring the_5cr sky n oise o f each image in the same 
fashion as lTeplitz et al.l (l2013lb the sky rms is measured 
in 1000 semi-random empty 51 x 51 pixel boxes, mul- 
tiplied by the correlat ion ratio for each mosaic from 
iFruchter fc Hook (2002;), and normalized to an aperture 
of 0.2"radius. All the data are drizzled to a pixel scale of 
0.03 arcsec pixel -1 , and therefore the correlated noise is 
highest in the NIR (largest pixels; 0.128 arcsec pixel -1 ), 
followed by the optical (0.05 arcsec pixel -1 ) and then 
the NUV (0.0396 arcsec pixel -1 ). The aperture choice 
of 0.2"radius is optimized for relatively compact galaxies 
in the optical and NUV, and slightly overestimates the 
depths of the NIR data. 

2.2. NUV Data 

The UVUDF observations were obtained in 2012 us¬ 
ing WFC3/UVIS in three filters, two different ob¬ 
serving mode s, and three orientations as described in 
iTenlitz et al.l ( 2013 b The first half of the data (referred 
to as Epochs 1 and 2) was obtained in 2 x 2 binned mode, 
while the second half (Epoch 3) was obtained in the un¬ 
binned mode, with the addition of “post-flash” to add 
internal background light to the image to mitigate c harg e 
transfer efficiency (CTE) degradation (see Sect ion 12.31) . 
While these data are described in detail in lTeplitz et all 
( 2013 b here we present improvements to the process¬ 
ing of the NUV data that significantly improve the final 
image mosaics. This includes a CTE correction, new as¬ 
trometric alignment, and new dark calibrations. In ad¬ 
dition to the improved darks and other improvements, 
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Table 1 

UDF Imaging Summary 


Instrument / 
Camera 

Filter 

Effect ive a 
Wavelength 

(A) 

Zero 

Point 

(ABMAG) 

Galactic 

Extinction 

(ABMAG) 

Number 
of orbits 

Exposure 

Time 

(B) 

5cr b 

Depth 

(ABMAG) 

Area 

(arcmin 2 ) 

PSF 

FWHM 

(arcsec) 

WFC3/UVIS 

F225W 

2359 

24.0403 

0.058 

16 

44,072 

27.8 

7.3 

0.11 

WFC3/UVIS 

F275W 

2704 

24.1305 

0.048 

16 

41,978 

27.8 

7.3 

0.10 

WFC3/UVIS 

F336W 

3355 

24.6682 

0.040 

14 

37,646 

28.3 

7.3 

0.09 

ACS/WFC 

F435W 

4317 

25.673 

0.033 

56 

134,880 

29.2 

12.8® 

0.10 

ACS/WFC 

F606W 

5918 

26.486 

0.023 

56 

135,320 

29.6 

12.8 e 

0.10 

ACS/WFC 

F775W 

7693 

25.654 

0.016 

144 

347,110 

29.5 

12.8 e 

0.09 

ACS/WFC 

F850LP 

9055 

24.862 

0.012 

144 

346,620 

28.9 

12.8® 

0.10 

WFC3/IR C 

F105W 

10550 

26.2687 

0.0079 

100 

265,459 

30.1 

4.6 

0.19 

d 





4 

8,100 

28.7 

8.2 f 


WFC3/IR C 

F125W 

12486 

26.2303 

0.0069 

39 

226,586 

29.7 

4.6 

0.19 

d 





5 

9,100 

28.6 

8.2 f 


WFC3/IR C 

F140W 

13923 

26.4524 

0.0048 

30 

82,676 

29.8 

4.6 

0.20 

WFC3/IR 0 

F160W 

15370 

25.9463 

0.0037 

84 

352,674 

29.9 

4.6 

0.20 

d 





5 

9,200 

28.2 

8.2 f 



Note. — UDF image mosaic information, including filters, zero points, extinction, exposure times, sensitivities, and area 
coverage for each of the eleven bandpasses. 

a Effective wavelength as calculated in ITokunaga &; Vaccal <120051 1. also known as the ‘pivot wavelength’. 

b Limiting 5cr depth base on sky noise in empty regions of the dete ct or in an aperture of 0.2 /; radius. __ 

c Deep IR data from UDF09 and UDF12 surveys (IQesch et al. 2010c b: Bou wens et al.I 2011b; Koeke moer et aO2 013: Elli s et aLl 
[20131) 

d Shallow IR data from CANDELS dGrogin et al.|[201U ; I Koekemoer et aLl 120 111) 
e Source catalog is trimmed to central 11.4 arcmin 2 , see Section EE3). 
f Source catalog trimming results in this covering 6.8 arcmin 2 , see Section im 


we applied a backgroun d subtract ion to each science ex¬ 
posure as described in iTenlitz et al.l (1201311 . to ensure 
that no leftover gradient is present in the final image mo¬ 
saics. These data have been released on the Barbara A. 
Mikulski Archive for Space Telescopes (MAST) website 
at http://archive.stsci.edu/prepds/uvudf/. They 
cover - 60% of the UDF FOV. 


2.3. Charge Transfer Efficiency 

The CTE of the WFC3/UVIS detector has signifi¬ 
cantly degraded over time, as radiation causes perma¬ 
nent damag e of the charge tra nsfer device (CCD) lat¬ 
tices (lMacKentv fe Smithll2012f) . This damage degrades 
the ability of electrons to transfer from one pixel to an¬ 
other, temporarily trapping electrons during the read¬ 
out. When uncorrected, the electrons are smeared out 
in the readout direction, appearing as trails in the im¬ 
ages. This affects the photometry and m easu red mor¬ 
phology of the ob je cts in the images (|Rhodes et al.ll20iot 
iMassev et alJl2010lh CTE degradation is most severe for 
low-background imaging of faint sources, such as NUV 
imaging and calibration dark frames, where faint sources 
or ho t pixels can be lost completely (lAnderson et al.1 
l2012fl . The effects of CTE degradation on the UVUDF 
data pro ducts w ithout any corrections are described in 
detail in ITenlitz et all (12013th who find that these is¬ 
sues are somewhat mitigated by the use of “post-flash”, 
thereby reducing the smearing effect, and avoiding the 
loss of faint sources and hot pixels. 

For this paper, only the post-flashed data are used, 
which significantly reduces the effects of CTE degra¬ 
dation. The post-flashed data are from “Epoch 3” 
and include ~13 electrons per pixel background, (see 
iTenlitz et al.1 I2013T ) . In addition, we apply a pixel- 


based CTE correctiorEB to the raw data based on 
empirical modeling of hot pixels in dark exposures 
(|Anderson fe Bedinll20 1 01 fMassev et al.ll20ldh . This cor¬ 
rection not only corrects the photometry, but also re¬ 
stores the morphology of sources. Currently, this soft¬ 
ware is only available for unbinne d data, a nd thus_the 
Epoch 1 and 2 data described in ITenlitz et al.l ( 2013 ) 
are not included in our analysis. While this reduces the 
overall depth of the NUV data, the CTE correction is 
critical to obtaining accurate photometric and morpho¬ 
logical measurements. CTE corrections are unnecessary 
for the optical and NIR data, as the ACS/WFC optical 
data were obtained early in the lifetime of the ACS/WFC 
and the HgCdTe NIR detector onboard WFC3/IR is not 
affected by CTE degradation. 

2.4. WFC3/UVIS Dark Calibrations 

Dark calibrations are especially important for NUV 
data because the dark current level in each exposure is 
high relative to the low sky background. In addition, 
regular calibration dark data can be used to identify hot 
pixels, which vary significantly over time. ITenlitz et al.1 
(120131) show that the darks currently provided by STScI 
are insufficient for data with low background levels after 
the CTE degradation of WFC3/UVIS. In this paper we 
improve the dark calibrations even further, as discussed 
below. 

While the STScI superdarks were mostly sufficient for 
early data obtained soon after the installation of WFC3, 
subsequent changes in the characteristics of the detector 
(such as CTE degradation) increasingly affected the sci¬ 
ence data. There are three major areas that the STScI 
processed superdarks are insufficient for use in the UDF 
program. First, when the darks produced by STScI are 
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used to calibrate the data, more than half of the hot pix¬ 
els are missed. Second, an observed background gradient 
is left unrenroved, and finally, a blotchy pattern is left in 
the images. We discuss these three effects in Appendix 
A. They are present even when the post-flashed raw dark 
data are corrected with the pixel-based CTE correction. 

In Appendix A we present a new dark processing 
methodolog y w hich impro ves o n the previous work by 
STScI and iTenlitz et al.l (12013T) . This new methodol¬ 
ogy is being used by a large number of other HST pro¬ 
grams (e.g., PI: Faber & Ferguson 12444, 12445, PI: 
Siana, 12266, 12201, 13389, PI: Guo 13309, PI: Malkan, 
12283, 12568, 12902, 13352, 13517), and we describe the 
methodology in detail here. A similar strategy is being 
developed at STScI, and will be implemented in their 
calibration pipeline in the future. 

2.5. Astrometric Alignment 

Here we discuss several sources of astrometric uncer¬ 
tainties in the original data, as well as our approaches 
to mitigating these and aligning all the images to a com¬ 
mon reference grid. The observations were all obtained 
in a non-integer pixel-offset dither pattern, aimed at en¬ 
suring that the point spread function (PSF) was ade¬ 
quately sampled in the final mosaics. Due to the geo¬ 
metrical distortion of the detector, shifts that are too 
large correspond to a substantially different number of 
pixels along the edge than at the center where the shifts 
are defined, which would cause the pixel subsampling 
phase to change across the detector. Therefore the shifts 
are kept small enough to retain the intended sub-pixel 
subsampling across the detector. Each small angle ma¬ 
neuver introduces a slight uncertainty in positioning (of 
the order of about 1 to 2 milliarcseconds). In addition, 
an optical offset is introduced when a different filter is 
inserted into the optical path. Moreover, during each 
orbit the spacecraft undergoes thermal expansion and 
contraction (“breathing”) due to changes in solar illumi¬ 
nation, which lead to changes in the optical path length 
to the detectors, hence resulting in slight scale changes 
from one exposure to the next. Finally, guide star reac¬ 
quisition uncertainties can lead to errors in position as 
well as small rotation uncertainties, while a full acqui¬ 
sition of a new guide star has astrometric uncertainties 
of ~ 0"3 — 0'/5 (reflecting the absolute astrometric uncer¬ 
tainties in the Guide Star Catalog 2). 

We make use of the source positions measured in 
the F435W ACS mosaics of the UDF ((Beckwith et al.1 
1200611 as our absolute astrometric reference frame. In 
this way, we minimize astrometric scatter resulting from 
color terms across sources. Initially we aligned the 
F336W to this frame, which then enabled us to boot¬ 
strap the F275W data, and subsequently the F225W 
data. The astrometric accuracy also depends on the 
accuracy of distortion calibration for WFC3; for this 
program, we found significantly improved results when 
using an updated distortion model made publicly avail¬ 
able by the WFC3 team which was implemented in 
the archive in 2013. The initial alignment was ac¬ 
complished with drizzlepac/tweakreg using catalog 
matching, which provides measurements of rotations 
as well as removing the bulk of the shifts that are 
present in the data. Further details of this tech¬ 
nique are presented in iKoekemoer et al.1 (1 20 1 ll) , updated 


with new distortion models that are available within 
drizzlepac/astrodrizzle. Once all the rotations and 
the majority of the shift information is solved for, the im¬ 
ages were then passed through astrometric refinement us¬ 
ing cross-correlation, which provides additional improve¬ 
ments, resulting in overall astrometric accuracies of ~2 
milliarcseconds in the mean shift positions of all the ex¬ 
posures relative to one another, which is the best possible 
level that is achievable given the sparse number of sources 
and their faint flux at these wavelengths, particularly in 
F225W. 

2.6. Extended UDF Near-IR Mosaics 

The UDF 12 program (|Koekemoer et all 120131 
lEllis et, all 1 201 31 ) originally released mosaics of the NIR 
UDF at the 60 mas scale, covering the extremely deep 
pointing of the NI R UDF from both th e UDF 12 and the 
UDF0 9 programs (lOesch et al.l 1 201 Od lbl: iBouwens et al.l 
1201 lbf) . The initial release of these mosaics do not cover 
the whole field of view of the original UDF or the new 
NUV observations (see Figure [[J. In addition, given the 
higher resolution and smaller plate scale of the NUV 
data, the NUV data mosaics are released at the 30 
mas scale to MAST, enabling impro ved morphological 
measurements in the UDF (e.g. JBond et al.ll2014f) . 

To facilitate photometry at the same plate scale as the 
NUV and optical data, we make use of new 30 mas scale 
NIR mosaics that include data from the surrounding 
CAN DELS survey (|Grogin et alJ[2CTTTt IKoekemoer et all 
1 201 ll ) as well as the full WFC3/IR dataset on the UDF 
footprint (IKoekemoer et al.l 120131 : lEllis et all 1201311 . In 
this way, the entire original optical UDF is covered by 
NIR data, although at a shallower depth in the regions 
outside of the UDF09/UDF12 NIR pointing. The depths 
and covered areas are listed in Table [1] The area of the 
shallower NIR data covers almost twice the area of the 
deep NIR pointing, a significant contribution to the UDF 
dataset. Moreover, while the data are shallower than the 
deep pointing, their depth is similar to our NUV data 
depths, and significantly improve the photometric red- 
shifts we describe below. 

2.7. Point Spread Functions 

In order to measure PSF-corrected aperture-matched 
photometry as described below, PSFs of the different 
bandpasses are required. PSFs for each HST camera 
are created in slightly different fashions, due to varying 
constraints of the data. 

The WFC3/NIR data PSFs are created in a very 
similar manner to the PSFs created by CANDELS 
(Ivan der Wei et alll2012H . Specifically, the NIR PSFs are 
hybrid PSFs, with the central pixels and the outer wings 
based on a model PSF, and the middle portion consisting 
of a stack of stars in the UDF. This hybrid approach is 
necessary since the large wings of the NIR PSFs contain 
approximately ~ 10% of the flux beyond 1"(WFC3 in¬ 
strument handbook), which is difficult to include solely 
by stacking the limited number of stars in the UDF. In 
addition, the NIR models were found to not accurately 
reproduce the PSF at intermediate radii compared to 
stars in CANDELS, there by justifying a hybrid PS F over 
using a pure model PSF (Ivan der Wei et a ill 20 121) . Note 
that the CANDELS PSFs can not be used directly, be- 
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cause we used a drizzle scale of 30 mas rather than 60 
mas. 

Th e NIR, P SF m odel is created with the TinyTim pack¬ 
age (iKristl 11995f t for the center of the WFC3 detec¬ 
tor, sub-sampled to align the planted PSFs, re-sampled 
to the WFC pixel scale, distortion corrected, and then 
combined with the same dither pattern and drizzle 
parameters as was used to produce the imaging mo 
saics. The resultant mo dels are similar to those from 
Ivan der Wei et all ( 2012 ). except at a 30 mas scale. The 
NIR stack of stars is created from 7 unsaturated stars 
in the UDF deep NIR region (Figure [lj se lected from 
the p ublished catalog of stars in the UDF dPirzkal et all 
l2005f) such that they are not contaminated in the wings 
by nearby galaxies. In addition, one star is excluded 
as it has a slightly extended profile and it is likely an 
AGN. The stars are registered to their subpixel centers, 
normalized by the peak value, and then coadded via a 
median. The final hybrid PSF is a combination of the 
two, composed of the PSF model from a radius of 0-5 
pixels (0-0.15"), of the stack of stars from 5 pixels to 35 
pixels (0.15-1.05"), and of the PSF model from 35-166 
pixels (1.05-5.0"). 

The WFC/ACS optical PSFs are median stacks of 
stars, constructed from 15 stars for the F435W and 
F606W bands, and 8 stars for F775W and F850LP bands 
selected from the pu blished catalog of stars in the UDF 
llPirzkal et all 12005 (1. The stars are registered to their 
subpixel centers, normalized by the peak value, and then 
coadded via a median. The optical bandpasses have a 
smaller PSF, with less flux in the wings and a narrower 
central peak, so we did not need to extend their PSFs to 
larger radii via a PSF model. 

The WFC3/UVIS PSFs are created with a hybrid PSF 
in a fashion similar to the WFC3/IR PSFs. The PSF 
FWHM of the NUV data is similar to the optical data, 
with < 5% of the PSF flux beyond 1"(WFC3 instrument 
handbook). However, due to the low number of stars 
in the NUV (2-3), the wings of the PSFs can not be 
recreated from a stack of stars. In fact, the best PSF in 
the NUV is generated by a single bright star. Coadding 
this higher SN star with the 1 or 2 low SN stars in the 
data does not improve the SN of the final PSF star in the 
NUV. Therefore, the best PSF is generated by combining 
a single high SN star with a PSF model. 

The hybrid NUV PSFs are created in a similar fashion 
as the NIR PSFs described above, with a few differences. 
The model PSF is only used in the wings, as the central 
part of the PSF is sufficiently constrained by the single 
star due to the higher resolution and smaller plate scale. 
At the same time, the model is incorporated at a smaller 
radius due to lower SN in the wings. In order to prevent 
discontinuities in the resultant PSF, the PSF model and 
the star are added together, weighted by a Gaussian with 
a full width of 15 pixels (0.45"). 

The FWHM values for all the PSFs are listed in Ta¬ 
ble [T] For this measurement, a symmetric Gaussian is 
fit to the final PSFs, although that is not a good fit 
to these non-Gaussian PSFs. However, the Gaussian 
fit only provides an estimate of the size of the FWHM 
for the table. The NUV PSF FWHM’s measured here 
are smaller than those measured bv lTeplitz et all (|2013T ) 
in the same d ata, which is expected because the data 
in iTeplitz et all (|2013fl were not yet corrected for CTE 


degradation, which artificially made the PSFs more ex¬ 
tended. The resulting PSF FWHM’s are similar to those 
measured bv iWindhorst et ahl (1201 111 , although slightly 
larger. While PSFs are created for all the images, only 
the F775W and the four NIR PSFs are used in this paper, 
as described in Section l3Jl 

3. PHOTOMETRY AND CATALOG CREATION 
3.1. Aperture-Matched PSF Corrected Photometry 


Aperture-matched PSF-corrected photometry is es¬ 
sential for robust photometry of galaxies using im- 


age data with varying PSFs (e.g.. Benitez et al. 

1999 

Vanzella et al.: j 2001l; 

Coe et al. 20061; Laidlcr et al. 

2007; 

de Santis et al. 2007 

: Finkelstein et al. 20121). As the 


PSF FWHM increases, the accuracy of galaxy photom¬ 
etry is affected. If sufficiently large apertures are used, 
then the photometry of bright galaxies are not signifi¬ 
cantly affected by this effect. However, for faint galaxies, 
a significant fraction of the flux is under the level of the 
noise, resulting in large differences in photometry regard¬ 
less of the aperture size. Multiple methods are possible 
to correct the photometry for the PSF variation, with 
the three most prominent ones discussed here. 

First, if the PSF FWHM’s are very different (by a fac¬ 
tor of a few or more), then a template-fitting method is 
typically used to match a convolved galaxy ‘model’ to th e 
measured flux (lLaidler et al.ll2007j : Ide Santis et al.l[2007h . 
Since our analysis is limited to data with similar PSF 
FWHM’s (within a factor of two), a simpler approach can 
be used. A second option is to convolve all the images 
to the PSF of the largest PSF image (e.g.. lGrazian et ahl 
120061 : [Finkelst ein et al.l [20l2h . and measure photometry 
on these PSF matched images. The downside of this ap¬ 
proach is that all measurements are then performed on 
convolved data rather than the original data, at the worst 
resolution of the images, and any systematic uncertain¬ 
ties in the PSFs will affect the photometry in all bands. 
This is particularly a problem for compact faint galaxies 
and low SN observations, for which the SN is signifi¬ 
cantly degraded once they are convolved by a larger PSF 
because the convolution includes pixels that are mostly 
noise. 

The third option, and the one used here, is to mea¬ 
sure the photometry in the original higher-resolution 
data, and then apply a PSF correction on the NIR mea¬ 
surements, which have larger PSF FWHMs, as done in 
iCoe et all 1)200611 . Since the PSFs of the NUV and opti¬ 
cal data are similar, PSF corrections between these band- 
passes could introduce systematic uncertainties based on 
the quality of the PSFs, which could be as large as the 
corrections themselves. More importantly, this method 
enables the use of a separate detection image for the 
NUV data as described in Section 13.61 which improves 
the SN in the NUV data by using smaller apertures, and 
afterwards applying an aperture correction. Since the 
emphasis of this paper is the inclusion of the NUV data, 
the third method is applied in this paper. 

3.2. Overview of Methodology 

We first summarize the general methodology, and 
then provide details in the following sub-sections. 
We use ColorPro to measure photometry in the im¬ 
ages ( Coe et aT] l2006t l. which is a wrapper for run- 
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Table 2 

Source Extractor Parameters for Different Runs 


Parameter 

Deep 

Shallow 

Deep Deblend a 

Shallow Deblend 3, 

NUV 

NUV Deblend a 

detect.thresh 

i.i 

3.5 

1.1 

3.5 

1.0 

1.0 

analysis.thresh 

i.i 

3.5 

1.1 

3.5 

1.0 

1.0 

deblend_nthresh 

32 

32 

8 

8 

32 

8 

deblend_mincont 

0.01 

0.01 

0.3 

0.3 

0.01 

0.3 

detect_minarea 

9 

9 

9 

9 

9 

9 

back_size 

128 

128 

128 

128 

128 

128 

back_filtersize 

5 

5 

5 

5 

5 

5 

back_photo_thick 

26 

26 

26 

26 

26 

26 


Note. — Table of SExtractor parameters used for the different runs to create the catalog. 
a Deblend represents the SExtractor runs with a low deblending threshold. 


ning Source Extractor (hereafter SExtractor) (v2.5.0; 
iBertin fe Arnoutsl 1 1 9961) on multiple images. Each run 
of ColorPro is analogous to running SExtractor in 
dual-image mode eleven times, once for each filter, and 
also applies aperture and PSF corrections (see Section 
13.31) . All photometry is therefore determined by running 
SExtractor in dual-image mode on both a detection im¬ 
age and each individual image. The detection image is 
created from multiple images obtained with different fil¬ 
ters, to maximize its depth and the robustness of the 
aperture sizes (see Section l3~fl) . 

A single detection threshold that optimizes the aper¬ 
tures of most of the objects in the images results in aper¬ 
tures that are poorly defined for both the bright sources, 
and fainter sources near them. The solution is to deter¬ 
mine apertures with two different detection thresholds. 
In addition, for each detection threshold, we also measure 
the photometry using two different deblending thresh¬ 
olds, as no single threshold sufficiently deblends some 
objects, without artificially splitting up single objects. 
This results in 4 separate measurements for each galaxy, 
with differing object definitions and aperture sizes, which 
are combined into a single photo metric catalog and seg¬ 
mentation map (see Section 13.51) . 

To optimize the SN in the NUV data, smaller apertures 
are needed, and are determined by running SExtractor 
in dual image mode with the F435W image, instead of 
the general detection image, for the three NUV images. 
Two separate deblending thresholds are required, and 
the resulting NUV photometry is merged with the full 
11-band catalog, including an aperture correction for the 
different aperture sizes (see Section 13.61) . 

3.3. Aperture and PSF corrections 

The PSF correction is determined by degrading the 
F775W image to the PSF of each of the 4 NIR images us¬ 
ing the IRAF task psfmatch. The F775W image is used 
as the high-resolution image, as it is the reddest high- 
resolution image with a wel l-behaved PSF (for a discus¬ 
sion of F850LP, see lCoe et al.ll2006lf . For each image, a 
PSF Kernel is created, which if convolved by the F775W 
image, matches the PSF of the corresponding NIR image. 
The psfmatch threshold parameter, which sets the low 
frequency cutoff in the Kernel creation, is set to 0.14 to 
minimize fringing in the resulting convolved image. 

ColorPro measures the colors of galaxies based on 
their isophotal fluxes, which have been shown to pro¬ 
duce robust col ors, outperform ing circular or large Kron 
apertures ijBenftez et al.1 120041) . It then applies an aper¬ 


ture correction from the ratio of the isophotal flux and 
the total flux measured via the SExtractor mag_auto 
Kron l|Kronl H980l) aperture flux. In addition, for the 
small number of sources in which an aperture and PSF 
correction cannot be calculated, ColorPro calculates it 
based on other sources with similar sizes (and thereby 
also magnitude). 

3.4. Multi-band Detection Image 

The detection image determines the object definitions 
and aperture sizes, and thus a very sensitive image is 
desired. It is also important for the apertures to be suf¬ 
ficiently large to contain the majority of the light in the 
worst resolution image, yet still be sensitive to smaller 
faint objects that may not be detected in the images with 
larger PSF FWHMs. A well-matched aperture minimizes 
the effects of any internal color gradients of the galaxies. 
Therefore, we make the detection image by averaging 
the 4 optical images with the 4 NIR images together, 
weighted by the inverse variance of each image on a pixel- 
by-pixel basis. In this fashion, the varied depth of the 
NIR image is taken into account in the region covered by 
the shallow NIR data, and the resultant RMS image is 
outputted for use with SExtractor. This detection im¬ 
age optimizes source detection at both resolutions as we 
do not lose the faint compact sources present in the op¬ 
tical data but not in the NIR data. Although averaging 
two different resolutions together creates slight ’donut’ 
images, this is sufficiently minor that it is not visible in 
the image, and SExtractor is not sensitive to this when 
defining the segmentation map as the detection image is 
convolved by a Gaussian in the SExtractor run. The 
alternative method of first smoothing all the data to the 
resolution of the largest PSF would reduce the SN of 
faint compact sources, losing these potentially interest¬ 
ing sources. This detection image methodol ogy is s imilar 
to that used by the CLASH team dPostman et alJ[2012l : 
ICoe et al.ll2013lf . 

3.5. Source Aperture Definitions 

While there are many options available to optimize 
SExtractor to obtain the best object definitions pos¬ 
sible, no single set of detection thresholds and deblend¬ 
ing parameters perfectly detects bright, faint, large, and 
small galaxies simultaneously. Either a lot of galaxies 
are blended together, or the larger galaxies are split 
into multiple segments. This is an intrinsic limita¬ 
tion with the current automation methods, although for 
the vast majority of sources, source definitions are ade- 


















quate. To solve these issues for the problematic sources, 
SExtractor is run four times on each image, with each 
run optimized differently, and then merged together into 
a final catalog as described below. Two are optimized 
via the detection threshold, and two include different de¬ 
blending thresho lds. This is simil ar to the ’hot’ and ’cold’ 
mode method bv lBarden et al.l (I2012D . with differences in 
the methodology being due to independent development. 
The different parameters are also summarized in Table 

m 

The first iteration is a ‘deep’ SExtractor run, with 
the detect_thresh and analysisjthresh. parameters 
set to l.lir, which determines the minimum deviation 
above the background RMS (including correlation cor¬ 
rections for the drizzle). A minimum of 9 contiguous 
pixels above this threshold are required to define a source 
as set by detectjninarea. The background sky for de¬ 
tection is computed by setting back_size to 128 and 
back_f iltersize to 5, and the local sky for photome¬ 
try is determined with a background annulus of 26 pixels 
(back_photo_thick). These parameters yield reasonable 
performance in detecting faint sources while minimizing 
spurious detections. 

With the above threshold, however, the source defi¬ 
nitions near bright targets are poorly defined for both 
the bright sources, and fainter sources near them. The 
solution is to run SExtractor with a higher detection 
threshold, with detect_thresh set to 3.5<r. This sec¬ 
ond ‘shallow’ run detects significantly fewer sources, al¬ 
though the bright source detections are well defined. The 
two catalogs are then combined by replacing the bright 
sources and their neighbors in the first run with the shal¬ 
low run measurements. Sources brighter than 22nd mag¬ 
nitude are replaced, which includes 55 bright sources and 
1014 fainter neighbors, which fall within 250 pixels of the 
brighter sources. Faint sources which are not detected in 
the shallow run within 250 pixels are also included, to 
avoid losing any sources. 

Both of these SExtractor runs have difficulty deblend¬ 
ing sources properly; either the deblending thresholds are 
set too low, and multiple objects are blended together, 
or the thresholds are set too high, and single galaxies are 
deblended into multiple sources based on knots or spiral 
structure in the galaxies. The solution implemented here 
is to run SExtractor twice; first with‘normal’ deblend¬ 
ing thresholds optimized for the vast majority of galaxies, 
with deblend_nthresh set to 32 and deblend_mincont 
to 0.01, and second with ‘low’ deblending thresholds op¬ 
timized for larger galaxies, with deblend_nthresh set to 
8 and deblendjnincont to 0.3. The threshold depen¬ 
dent source definitions determined in the first two runs 
are manually checked for galaxies that should not be de¬ 
blended, and replaced by these 3rd and 4th runs. This 
only affects 90 galaxies in the final catalog, although 
a larger number of galaxies are deblended and should 
not be. This is due to the conservative approach used, 
which only combined galaxies that are obviously incor¬ 
rectly deblended to avoid combining multiple individual 
galaxies together. A small number galaxies can not be 
combined together with this technique, even when clearly 
they should be, as other nearby galaxies would be in¬ 
correctly deblended in those cases. Hence, some galaxy 
blending and deblending issues will be in the final cata¬ 
log. The resulting catalogs are all merged into a single 


catalog, with a single segmentation map. This is then 
fed back into ColorPro to create a final catalog ready 
for calculating photometric redsliifts. 

3.6. F435W Apertures For NUV Photometry 

One of the primary goals of this pa per is to inclu de th e 
new NUV observations of the UDF (ITeplitz et al.ll20l3l i 
in the galaxy photometry. While the gala xies in t he dat a 
have similar rest-frame optical sizes ([Bond et al.1120141 ). 
the images are somewhat shallower than the rest of the 
UDF data (see Tabled, and the galaxies in the NUV are 
fainter due to their spectral energy distribution (SED). 
While the apertures determined in Section[X5]could work 
for the NUV data, the SN is not optimized when using 
very large apertures on small faint objects. In addition, 
the use of smaller apertures helps avoid remaining cali¬ 
bration issues observed as a blotchy pattern in the dark 
calibrations which results in an increased flux and uncer¬ 
tainty over very large areas (see Section IAT sI) . Therefore, 
smaller apertures are used for the NUV data determined 
from the F435W mosaic of the UDF, enabling more ro¬ 
bust and higher SN measurements. 

SExtractor is run in dual image mode with the F435W 
mosaic as the detection image for all three NUV mosaics 
and the F435W image. Since sources in the F435W im¬ 
age are fainter than in the redder bands, a single thresh¬ 
old selected such that no sources are detected in the 
negative image is sufficient, with detect_thresh=1.0<7 
and detect_minarea=9. By requiring no sources are de¬ 
tected in the negative image with the same SExtractor 
parameters, we ensure that the parameters are optimized 
since the negative image has the same noise properties 
as the normal image, but with no real sources. In this 
way there will only be almost no false sources detected 
by SExtractor when run on the original i mages, except 
by potential image artifacts (jRafelski et al.ll2009f ). 

Galaxies in the F435W band are clumpier and have 
a hi g her l evel of apparent structure in the images 
ijGuo et al.l 120151) . making two runs with different de¬ 
blending thresholds more important than in the multi¬ 
band detection image. The same two deblending thresh¬ 
olds are used as in Section [3751 with a normal and low de- 
blending threshold. Sources detected in the low deblend¬ 
ing threshold case are cross-checked with the merged 
multi-band detection segmentation map from Section 
13.51 If the segmentation map contains more than a single 
segment in the area covered by the smaller F435W seg¬ 
ment area, then the normal deblending threshold is used 
as otherwise multiple galaxies are likely to be blended to¬ 
gether. Otherwise, the low deblending threshold is used, 
to maximize the combination of clumps and other struc¬ 
ture in galaxies into single galaxy detections. In this way, 
the F435W source definitions are similar to the source 
definitions of the multi-band detection definitions. 

The new F435W detected sources are then matched to 
the multi-band detected sources using the merged seg¬ 
mentation maps described above. For each source in the 
multi-band catalog, all the pixels in the aperture defined 
by the segmentation map are checked in the F435W seg¬ 
mentation map. If only a single source is found within 
that aperture, then the source is matched. If two sources 
fall within the aperture and the second object covers 40% 
or more of the area of the source, then the two sources are 
added together and errors propagated. Otherwise, only 
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the main source is matched to the multi-band catalog. 
If no source is detected in the F435W band within the 
detected area, it is marked as undetected in the catalog. 
Since the NUV data are significantly shallower than the 
F435W band, all sources not detected in F435W are as¬ 
sumed not detected in the NUV bands, which will almost 
always be the case. This methodology works well because 
if there are multiple sources in the aperture, then they 
were detected with the low deblending threshold param¬ 
eters, and most galaxies are either one or two sources 
before being matched. In addition, if a source is close to 
the edge or the chip gap of the NUV data, it is marked 
as not covered. 

Since the F435W detected apertures are smaller than 
the multi-band apertures, an aperture correction is added 
to place it on the same magnitude scale as the larger 
aperture longer wavelength data. The aperture correc¬ 
tion is determined from the higher SN F435W image, 
consisting of the difference in the large aperture and 
small aperture magnitudes. As with all aperture cor¬ 
rections, we assume the color gradient from F435W to 
the NUV is small. 

4. PHOTOMETRIC REDSHIFTS 


Photometric redshifts are an accepted and ro¬ 
bust procedure to estimate redshifts of galaxies 


without spectroscopic 
Benftezl 12000: Coe et al. 

information (e.g., Koo 
20061 Ilbert et al, 2006, 

1985; 

2009; 

Hildebrandt et all 20101 

Dahlen et al. 2013(1. There are 


a large number of photometric redshift software packages 
available, a nd the differences betwe en th em are sum¬ 
mariz ed by iHildebrandt et all ([2010) and iDahlen et al.l 
(2013). In addition to differences in the methodologies 
of the code, the choice of templates and filters affect the 
final redshift estimates. For these reasons, we present 
and compare redshifts from two of the best-performing 
packages with different methodologies, priors and tem¬ 
plate sets. 

The Bayesian Photometric Redsh i ft (BPZ) algorit hm 
(iBenftezl 120001 : iBenftez et all 120041: iCoe et aD 1200611 is 
the primary r edshift software used d ue to its robust 
performance (Hildebrandt et all I2010H . familiarity of 
the code to the authors, and best performance as 
detailed in Section 14.21 The BPZ SED templates 
and priors have been significantly u pdated (Benitez 
et a l. in prep.) as described in ICoe et all fl 201 .'II I 
and IBenftez et al.l (12014H . Specifically, the model spec¬ 
tral energy distr ibutions (SEDs) are base d on those 
from PEGASE (|Fioc fc Rocca-Volmerangel I1997H but 
re-calibrated based on observed photomet ry and spec¬ 
trosc opic redshifts from FIREWORKS (iWuvts et al.1 
2008[). The templates are shown in Figure 10 of 
Benitez et al.l ( 2014 ). and in Figure QJ| here. These tem¬ 
plates include four elliptical galaxies, one Lenticular, 
two spirals, and four starbursts. The templates in¬ 
clude emission lines, and we interpolate between each 
pair of adjacent templates to create nine intermedi¬ 
ate templates, yielding 111 possibilities in BPZ. The 
Bayesian prior is based on luminosity functions ob- 


served in COSMOS 

Tlbert et al 120001. GOODS-MUSTO 

(Grazian et al. 2006 

Santini et al. 2009i), and the UDF 

(Coe et al.j 2006[). The software is also modified from the 


publicly released version with minor improvements, such 
as matching the output uncertainties to the observed 


scatter. 

In addition to the BPZ redshifts, we also determine 
a seco nd set of reds hifts using the EAZY software 
dBrammer et al.l l2008fl . another top performer by the 
IHildebrandt et al!T ( 2010fl study. This code not only pro¬ 
vides another independent meth odology including th e 
use of a template error function dBrammer et al.ll2008f ). 
but also uses a different set of galaxy templates in- 
cluding emission line s based on star formation rates 
(iBrammer et al.lf201 1), and the inclusion of a very dusty 
template (Brammer et al. in prep). 

The quality of the photometric redshift fit is described 
by two quantities, ODDS and reduced chi square. The 
ODDS parameter measures the spread in the probability 
distribution function, P(z ), representing the probability 
of a galaxy being at a specific redshift. The ODDS for 
BPZ is defined as the integrated P(z ) contained within 
2*0.03*(l + z). When P{z ) has multiple peaks, the value 
of ODDS is low and the resulting redshift uncertainties 
are large. If P{z) has only a single peak, then the ODDS 
describes the width of the distribution, with a maximum 
value of ODDS of 1.0. Typically 0DDS> 0.9 is considered 
to be a high-probability determination of the re dshift 
(Bc mtedl2000l : IBenftez et aLll2004L ICoe et aPl2006H . 

The second quantity is a reduced chi square goodness 
of fit. In addition to the regular chi square used in the 
calculations, BPZ also reports a modified reduced chi 
square, Xmod- The Xmod similar to a normal reduced 
chi square, but includes an additional uncertainty for 
the SED templates in addition to the uncertainty in the 
galaxy photometry (ICoe et al.ll2006f) . The added uncer¬ 
tainty for the SED was determined such that the resul¬ 
tant Xinod * s a more realistic measure of the quality of 
the fi t (for more discussion of Xmod> see Rafclski et al. 
[2009 ). For the EAZY algorithm, a normal reduced chi 
square, xt: is used. While strictly speaking a xt ~ 1 
would be considered a good fit, in practice a Xmod < 4 
was foun d to indicate a relat ively good fit with robust 
redshifts (jRafelski et al.ll2009 ). 

A comparison of Xmod f rom BPZ and xt from EAZY 
finds that they differ as a function of magnitude, with 



Figure 3. BPZ templates including four elliptical galaxies (Ell), 
one Lenticular (ESO), two spirals (Sbc and Scd), and four star- 
bursts (SB). Th ese templates ar e based on those from PEGASE 
JFioc &; Rocca-Volmerangd 119971) but re-calibrated based on ob¬ 
served photometry and spectroscopic redshifts from FIREWORKS 
(IWuvts et al.|[2008l) . 
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Table 3 

Spectroscopic Redshift Compilation 


Survey a 

Telescope* 5 

Instrument 0 

Number 

Quality 

Reference 

VVDS 

VLT 

VIMOS 

8 

95% 

Le Fevre et al. (2004) 

Szokoly 

VLT 

FORS1/FORS2 

4 

‘reliable’ (2 or 2+) 

Szokoly et al. (2004) 

K20 

VLT 

FORS1/FORS2 

27 

‘secure’ (1) 

Mienoli et al. (2005) 

GRAPES 

HST 

ACS 

6 

UV absorption lines 

Daddi et al. (2005) 

Vanzella GOODS 

VLT 

FORS2/VIMOS 

39 

A or B multiple lines 

Vanzella et al. (2005, 2006, 2008, 2009) 

Fopesso GOODS 

VLT 

VIMOS 

19 

A or B multiple lines 

PoDesso et al. (2009) 

Balestra GOODS 

VLT 

VIMOS 

18 

A or B multiple lines 

Balestra et al. (2010) 

GMASS 

VLT 

VIMOS 

29 

‘good’ (1) 

Kurk et al. (2013) 

3D-HST 

HST 

WFC3 

28 

3 or 4, multiple lines 

Morris et al. (2015) 


Note. — Table describing spectroscopic redshift compilation. Only robust redshifts are used, with multiple lines identifying 
the redshifts. 

a VVDS stands for VIMOS VLT Deep Survey, GRAPES for Grism ACS Program for Extragalactic Science, (IPirzkal et al.|[2004l) , 
GOODS for Great Observatories Origins Deep Survey, G MASS for Galaxy Mass Assembly ultra-deep Spectroscopic Survey, and 
3D-HST is a large HST program (|B rammer et al.l [20121 ). 
b Either Very Large Telescope (VLT) or Hubble Space Telescope (HST). 

c Either FOcal Reducer and low dispersion Spectrograph (FORS), Visible Multi-Object Spectrograph (VIMOS), Advanced 
Camera for Surveys (ACS), or Wide Field Camera 3 (WFC3). 


x1 being larger. For fainter galaxies (F606W magnitude 
> 27), they differ by a factor of ~ 2. However, as the 
brightness of the galaxies increases, the ratio increases as 
does the scatter. For bright galaxies (F606W magnitude 
< 25), the ratio increases to ~ 11. This is reasonable, 
given that xiLod was specifically designed to correct for 
that fact that galaxies with small photometric un cert ain¬ 
ties tend to have high y? even for good fits (|Coe et al.l 
l20f)3 IRafelski et a,Tl20F)9 ) , since systematic uncertainties 
for the SED templates are not taken into account. In or¬ 
der to have a similar number of targets selected as having 
good x 2 f° r both BPZ and EAZY in the investigations 
below, a xt < 10is used as an indicator of a good fit for 
EAZY redshifts in this paper. Restrictions to lower y 2 
and higher ODDS will improve the redshift quality further, 
yet at the cost of sample size. 

While naively one may expect that the inclusion of 
mid-IR data may improve the photometric redshifts, em¬ 
pirical tests show that many algorithms (including BPZ) 
actually perform worse when inc luding the mid-IR pho¬ 
tometry (Hildcbra ndt et aLll2010D . This was confirmed 
in tests on our data set as well, and is likely due to in¬ 
sufficient knowledge of the template SEDs in the mid- 
IR, due to degeneracies in the dust and PAH spectral 
features. In addition to such mid-IR data having sig¬ 
nificantly worse PSF FWHM’s than the data described 
here which complicates its use, and the data being signif¬ 
icantly less deep than the other data sets used here, there 
is also no improvement in the photometric redshifts with 
their inclusion. Therefore, such data are not included for 
the photometric redshift determinations. 

The catalog and photometric redshifts are not opti¬ 
mized for searches of z > 8 galaxies, due to choice made 
in the catalog creation to optimize intermediate redshift 
source selection. First, as described in Section 13.41 the 
detection image consists of 8 images including the optical 
data. Hence, faint galaxies appearing only in the reddest 
NIR bands may not be detected at all, as the NIR flux is 
averaged with that in bluer bands in which the galaxy’s 
Lyman break flux decrement causes the galaxy to drop 
out. Furthermore, as discussed in Secton[5l single band 
detections are excluded from the catalog to avoid spu¬ 
rious detections, and hence the highest redshift galax¬ 


ies that only appear in the F160W image are excluded. 
While we find no new z > 8 galaxies, we do recover can¬ 
didates out to z ~ 8.3 and NIR mag ~ 29 previously 
published in the UDF (lEllis et al.1 120131 iMcLure et all 
120131 : ISchenker et al.ll2013HOesch et al.l 1201 T 

4.1. Sample of Spectroscopic Redshifts 

In order to assess the quality of the photometric red¬ 
shifts, a large sample of high-quality spectroscopic red¬ 
shifts is needed. We compiled a total of 176 robust spec¬ 
troscopic redshifts from 9 sep arate su rveys, which almost 
doubles the compilation from IRafelski et all (120091) . Al¬ 
most half the resultant sample have more than one red¬ 
shift measurement, and three of the measurements are in¬ 
consistent. For those three, the best redshift is selected 
based on the quality of the original data. While the 
number of spectroscopic redshifts could be significantly 
increased by including lower-quality redshifts, doing so 
would introduce an uncertainty when photometric and 
spectroscopic redshift measurements do not agree, mak¬ 
ing it unclear which redshift is incorrect. 

The source of the spectroscopic redshifts are outlined 
in Table [3J along with the telescopes and instruments 
used, the number of sources from each survey, and the 
quality cut applied. In general, the quality cut requires 
redshifts to be identified with multiple lines and good SN. 
The 3D-HST redshifts used differ from the official 3D- 
HST grism redshifts, which include photometry in their 
redshift determinations (see Sectio n 14.3D , and are from 
and independent analysis bv iMorris et al.l (I2015D . None 
of the redshifts used include photometry to determine 
the redshift, so they are an independent check on the 
photometric redshifts. Of the 176 spectroscopic redshifts 
described above, 7 are stars and are not included in the 
investigations below. 

The 169 reliable spectroscopic redshifts of galaxies are 
listed in Table [I] The table includes the main catalog 
ID number, spectroscopic redshift, photometric redshift, 
ODDS, and Xmod- Although not included in this table, the 
photometric redshifts of the 7 stars are incorrect because 
a stellar template is not included in the pho tometri c fits, 
and are all identified in the star catalog (IPirzkal et al.l 
l2005f ). Figure Q] shows a histogram of the galaxy red- 











































11 


Table 4 

Reliable Spectroscopic Redshifts of Galaxies in the UDF 


ID 

RA 

DEC 

^spec 

2 BPZ a 

xLa b 

0DDS c 

d 

Z EAZY 

A+ 

0DDS f 

Reference 8 

373 

53.15437485 

-27.82148101 

1.14 

i 99+0.10 
± ‘^-0 .10 

0.51 

0.99 

1 11+ 012 

1.22 

1.00 

5 

534 

53.16170093 

-27.81925383 

0.67 

U.DD — o os 

0.10 

1.00 

u - ol -0.11 

0.56 

1.00 

6 

865 

53.17450843 

-27.81495550 

0.67 

fi 64+ 0 - 07 
u -°^-0.08 

0.06 

0.99 

0 Ri+°- 09 
u - ol -0.10 

0.84 

1.00 

1 

983 

53.14989262 

-27.81400031 

1.31 

1 90+O.IO 

1.^8_0 10 

0.22 

0.99 

1 24”*” 011 

1.44 

1.00 

6 

1035 

53.17634676 

-27.81475285 

2.44 

9 49 + 0.15 
Z ' 4Z -0.15 

0.68 

1.00 

1 QQ+ 013 
i.yy_i 94 

7.85 

0.93 

6 

1060 

53.15915906 

-27.81374038 

1.77 

1 74 + 0 - 11 
'^-0.13 

0.27 

0.99 

1 71 +0-08 

J " ' ± —0.08 

2.45 

1.00 

8 

1077 

53.16529350 

-27.81405511 

3.06 

0 09+O.I8 
°'°^-0.19 

0.73 

1.00 

3 04”*" 0-29 

°- u ^-0.17 

1.03 

1.00 

2 

1134 

53.16817307 

-27.81293086 

0.96 

1 no” 1-9 ' 99 

l.UU — o 09 

0.18 

0.99 

0 7Q+0.09 
u - /y - 0.10 

1.48 

1.00 

6 

1220 

53.17932560 

-27.81252367 

1.77 

1 qq+0.12 

1 . 8 cS_o 13 

0.89 

0.99 

1 74+ 0 06 
± ‘ * ^—0.06 

11.48 

1.00 

9 

1438 

53.14508947 

-27.80985318 

1.25 

1 2 fi + 0 ' 09 
1 ,ZD — 0.11 

0.01 

0.99 

1 9-1+O.IO 

± - z± -o.n 

0.93 

1.00 

9 


Note. — Table [4l is published in its entirety in a machine-readable form in the online version of the Astronomical Journal. A portion 
is shown here for guidance regarding its form and content. 

a Bayesian Photometric Redshift (BPZ) and uncertainty from 95% confidence interval. 

Modified reduced chi square fit, where the templates are given uncertainties. 
c Integrated P(z) contained within 2 * 0.03(1 + 2 bpz)- 
^ EAZY redshift and uncertainty from 95% confidence interval. 
e Reduced chi square fit. 

^ Integrated P(z) contained within 0.2(1 + zeazy)- 

8 (1) Le Fevre et al.l J2004D (2) ISzokoly et al. (2004) (3) Mignoli ct al. (2005) (41 D addi et all J2005D (5) iVanzella et al.l J2005I , [20061 , 
1 20081 , 20091) (6)[Popesso ct al. (2009) (7) Balestra ct al. (2010) (8) Kurk ct al. (2013) (9) Morris et al. ( 2 0151) 



Figure 4. Spectroscopic redshift distribution of 169 galaxies in 
the Hubble Ultra Deep Field. Most of the spectroscopic redshifts 
are at 2 < 2 . 


shifts, and shows that the vast majority of spectroscopic 
redshifts are at z < 2, with a peak at z ~ 1. This 
means that the photometric redshifts can be well tested 
at z ~ 1, but that redshifts at 2 > 2 will be less well 
vetted. Some improvement on the redshift distribution 
is obtained with the grisrn redshifts, although with their 
own caveats as discussed in Section 14.31 
Neither of the redshift codes is explicitly tuned or 
trained using these spectroscopic redshifts, although the 
templates used by BPZ are re-calibrated based on ob¬ 
served photometry and spec t roscopi c redshifts from the 
FIREWORKS catalog dWuvts et alJl2008l) . This catalog 
includes about half of the redshifts presented here from 
spectroscopic campaigns before 2008. This does not sig¬ 
nificantly affect our quality checks below, as there are 
over 6000 redshifts in the FIREWORKS catalog which 
are used in refining the templates. Moreover, the pho¬ 
tometry used would not include the new NUV or NIR 
data in the UDF, which are primarily responsible for the 
improvements in the redshifts derived here. Lastly, the 
improvements in redshift derivations with the addition of 


the NUV and NIR data discussed in Section are com¬ 
pared to redshifts determined with the same templates, 
and thus the comparisons are unbiased. 


4.2. Quality of Photometric Redshifts 

The quality of the photometric redshifts is evaluated 
by a comparison of the subset of galaxies with high qual¬ 
ity spectroscopic redshifts, such that it is extremely likely 
that any discrepancies are due to inaccuracies of the pho¬ 
tometric rather than the spectroscopic redshifts. Figures 
[5] and [6] show the very good quality of the photomet¬ 
ric redshifts in the UDF, with very few outliers and a 
tight relation around the one-to-one line. The compari¬ 
son plot is on a ‘pseudolog’ (log(l + z)) scale to spread 
out the data points at lower redshifts. No optimizations 
are performed for this comparison, so it is representative 
of the general quality of the photometric redshifts, al¬ 
though the galaxies with spectroscopic redshifts consist 
of the brighter galaxies in the UDF. The quality of the 
redshifts are quantified by the scatter in the difference be¬ 
tween the photometric and spectroscopic redshifts, and 
by the fraction of outliers far from the line. 


4.2.1. Scatter of Redshifts 

The scatter in the difference of the photometric and 
spectroscopic redshifts provides a metric to quantify the 
general accuracy of the photometric redshifts. To min¬ 
imize the impact of outliers on the measurement of 
the scatter, the no r malized m edian absolute deviation 
f(TMMin: lBrammer et alJl2008h is used to define the scat¬ 
ter: 


onmad = 1-48 x median 


A: — median(Az) 


1 + z 


spec 


(i) 


where A z = z spec — z p hot, -z sp ecis the spectroscopic red¬ 
shift, and z p hot is the photometric redshift. The onmad 
is equivalent to the standard deviation for a Gaussian 
distribution, but is less sensitive to outliers. There is a 
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Figure 5. Photometric versus spectroscopic redshifts for 169 
galaxies in the UDF on a ‘pseudolog’ scale, showing the high quality 
of the BPZ photometric redshifts. The blue circles represent red- 
shifts for galaxies that are covered by the NUV data, and the green 
squares are for redshifts without NUV data. All redshift determi¬ 
nations include NIR and optical data. Filled symbols represent 
photometric redshifts with good ODDS and chi square, ODDS >0.9 
and Xmod < 4 for BPZ and ODDS > 0.9 and xi < 10 for EAZY. The 
open symbols on the other hand are for redshifts that may have 
multiple peaks in P(z) or may have a low quality fit, ODDS < 0.9 or 
X^d > 4 for BPZ and ODDS < 0.9 or xj > 10 for EAZY. Dashed 
lines represent outlier cutoff, defined as \Az\ /(I + z spec ) < 0.15. 



Figure 6. This figure is the same as Figure |5] except for the 
EAZY photometric redshifts instead of the BPZ redshifts. 


slightly differ ent version of o nm a d also used in th e lit- 
erature fe.g.. Illbert et all 120091 : [Pahlen et al.ll2013fl . but 
the above version is more robust to outliers due to the 
use of the second median. 

The ctnmad is very small for our sample, significantly 
improving on previous photometric redshifts in the UDF, 
with onmad = 0.028 for BPZ and ctnmad = 0.030 
for EAZY. These numbers drop even further for BPZ 


when considering only the sample with NUV data, with 
onmad = 0.026, while there is an increase for EAZY, 
with onmad = 0.035 (see Section [4|. 

The las t full UDF photometric redshift catalog by 
iCoe et ah! (120061 1 has a onmad =0.068 using the same 
spectroscopic sample, and therefore the redshifts pre¬ 
sented here are a factor of two better than before, based 
on this metric. Another catal og o f photometri c red shifts 
in the UDF is presented bv iCameron et all (l201l!l . al¬ 
though it only includes the 1052 galaxies brighter than 
magFi60W < 27 that are in the smaller deep NIR region. 
Matching this sample to the spectroscopic redshifts de¬ 
fined above yields onmad =0.042 based on 93 redshifts. 

In addition to the UDF redshift catalogs, two recent 
catalogs of GOODS-South encompass the UDF, and have 
photometric redshifts for all 168 spectroscopic redshifts. 
The CANDELS team calculated redshifts by combining 
multiple redshifts from different software packages, us¬ 
ing 14 bandpasses including ground-based u-band data 
and Spitzer data, and have a onmad =0.022 (Santini 
et al. 2015, in prep). The 3D-HST team also calcu¬ 
lated photometric redshifts using EAZY, using up to 
40 bandpasses including significant numbers of ground- 
based data and Spit zer data, and have a onmad =0.013 
ijSkelton et al.ll2014ll . While these new catalogs have bet¬ 
ter onmad than our catalog, our results have about a fac¬ 
tor of two improvement in outlier fraction (see Section 
14. 2. 211 below. Our onmad also compares well to other 
larger studies of different fields with many bandpasses, 
such as the COSM OS redshifts (onmad = 0.007 — 0.033; 
Illbert et~aJ1l2009T) or the C ANDELS GOODS- south red¬ 
shifts (onmad = 0.03 — O.U lDahlen et al.ll2013l) . although 
our galaxy sample is fainter. 

4.2.2. Outlier Fractions 

Photometric redshifts can have “catastrophic” redshift 
errors, caused by template mismatches, when the wrong 
redshift is assigned to a galaxy, and the probability distri¬ 
bution function of the photo metric reds hift does not in¬ 
clude the true redsh i ft ( |EUis||1997|: |Ferna ndez-Soto et all 
119991 : iBenfted 120001 : IRafelski et all I2009H . Most of the 
time the probability function correctly includes both 
redshift possibilities, but if there is a single peak at 
the wrong redshift with an acceptable chi square value, 
then the galaxy is unambiguously assigned to the wrong 
redshift. The rate of these catastrophic errors in the 
UDF has signifi c antly decreased over previous work (e.g., 
ICoe et aD 120061 : IRafelski et al.ll2009H . but still exists for 
various reasons. 

Any individual galaxy redshift is only as good as its 
object definition, and bad object definitions will result 
in poor photometric redshifts. For instance, if multiple 
galaxies are included in a single aperture due to insuf¬ 
ficient deblending, then the resultant photometric red¬ 
shift will be incorrect. While every attempt is made 
to avoid this in the catalog, one such example is very 
evident in the comparison in the top left panel of Fig¬ 
ure 0 Object 10157 is composed of two overlapping 
galaxies at high and low redshifts, which results in an 
incorrect photometric redshift at z p hot ~ 2 in the Fig¬ 
ure. While the cause of the incorrect photometric red¬ 
shift is understood, it and others like it are still in¬ 
cluded in all subsequent analysis to avoid any selec¬ 
tion biases. A random selection of galaxies using pho- 
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tometric redshifts would likewise also include such ob¬ 
ject definitions. Therefore, including these galaxies in 
the comparisons below makes the outlier fractions con¬ 
sistent with what a user of the catalog will experience. 
For small samples of galaxies, users may wish to manu¬ 
ally inspect the segmentation map and the photometric 
redshift fits, which would yield a more robust sample. 
Galaxy cutouts and SED plots are available on the web 
at: http://www.rafelski.com/uvudf/catalogs.html 
for the BPZ redshifts. 

There are multiple methods to define an outlier in a 
comparison of spectroscopic and photometric redshifts. 
One method is to consider objects at > 5ctnmad as out¬ 
liers, since these r edshifts are significant ly more deviant 
than the scatter dBrammer et al.1 |2008f) . However, in 
that case the photometric redshifts with a larger scat¬ 
ter and the same number of significantly deviant red¬ 
shifts may have a lower outlier fraction than a dataset 
with a smaller scatter. An alternative method is compare 
the outliers t o an absolute deviation, such as |Az| /(I + 
z sr ,rr) > 0-15 (Ebert et al.1 20061 l2009t iHildebrandt et al.l 

120101 iBrammer et al.ll201l lDahlcn o t al.ll20 1 .'ill . although 

this could result in false outliers if the scatter is large. 
Given the onmad of our sample, these two methods are 
almost identical for our photometric redshifts, and the 
second method is used for the outlier fraction (OLF) dis¬ 
cussed below. 

There are a total of 4 outliers with |Az| /(I + z spec ) > 
0.15 for BPZ redshifts, and 10 outliers for the EAZY 
redshifts. These correspond to OLFs of 2.4% and 5.9% 
respectively. If limiting the spectroscopic sample to those 
with with good ODDS and y 2 as defined in Sectional then 
the outlier fractions change to 2.4% and 3.9% for BPZ 
and EAZY respectively. The OLF for EAZY drops to 
2.7% when considering only the sample with NUV data 
(three outliers), while the BPZ fraction stays the same at 
2.4% (although with one fewer outlier). The OLFs can 
be made even smaller if a more stringent cutoff for y 2 is 
used, although that also reduces the size of the sample. 

This catalog is a significant improvement over the pre- 
vious catalogs of the UDF, which have an OLF of 16.4% 
(|Coe et al.lI2006T) and 4.3% dCameron et al.l[2 01 ll) when 
determined with the new larger spe ctrosco pic s amp le de¬ 
fined above. Note that the OLF of lOoe et all ( 2006 ) de¬ 
creases to 10.6% when only including galaxies with good 
ODDS and y 2 . Our redshifts thereby provide a factor of 
three improvement in outlier rate than the previous full 
redshift catalo g of the UDF, and a factor of two better 
than those of ( Gameron et alJfeOll ). 

The catalog also compares well to the GOODS-South 
catalogs, which have outlier fractions of 3.0% and 5.3% 
for CANDELS and 3D-HST respectively when compar¬ 
ing to the same 168 spectroscopic redshifts. When com¬ 
paring to the smaller subsample that also have NUV 
data, both the CANDELS and 3D-HST outlier fractions 
stay about the same at 3.3% and 5.6%, showing no im¬ 
provement in this subsample of spectroscopic redshifts. 
This suggests that it is indeed the addition of the NUV 
data, and not sample selection, that is responsible for the 
improved outlier fractions in the NUV sample above. 

4.2.3. Combining Multiple Redshift Methods 

These outlier galaxies are not the same in the two pho¬ 
tometric redshift samples, and thus outliers can be fur¬ 


ther avoided by considering only galaxies with similar 
photometric redshifts between the two software pack¬ 
ages, defined as \Az p hot\ /(l + zbpz) < 0.15, where 
Az p hot = zbpz — ~eazy- The sample that consists of 
matching photometric redshifts yields an OLF of 0.63%, 
which consists of a single outlier. Restricting the sample 
to good ODDS and y 2 for both redshift codes keeps the 
single outlier, with an OLF of 0.68%. Requiring match¬ 
ing redshifts between codes reduces the sample size to 
93% and 86% for the two respectively. 

The matched photometric redshift sample does not sig¬ 
nificantly improve the scatter onmad, with onmad = 
0.027 for both BPZ and EAZY when restricting to 
matche d pho t omet ric redshifts. On the other hand, 
iDahlen et al.l (|2013f) found that combining multiple red¬ 
shifts results yielded the best photometric redshifts. Tak¬ 
ing an average photometric redshift from the two codes 
improves onmad for the sample with good ODDS and y 2 
for both methods, with onmad = 0.026 for the average 
redshift value. 

Since the improvements of combining both codes helps 
improve outlier fractions, and this is another option for 
the reader if they so choose. However, for simplicity the 
rest of the paper examines only a single redshift code. 
While both redshift codes provide similar quality red¬ 
shifts, the BPZ redshifts are utilized henceforth due to 
slightly better onmad and OLF (see also Section 14.3|) , 
and the authors familiarity with the code. The catalog 
in Section [5] presents both redshifts, and hence the reader 
can choose which redshifts to use. 

4.3. Comparison with Grism Redshifts 

The use of the WFC3 grisms on HST to d etermine red¬ 
shifts is increasing a t a rapid pace (e.g .. [Brammc r et al.1 
l2012t iGolbert et al.1120131 : lAtek et al.112014H . providing 
large numbers of galaxy redshifts. The UDF is cov¬ 
ered by the G141 WFC3 grism by both 3D-HST and 
CANDELS, and a special data release with redshifts 
(|Brammer et al.l[2012ll20T3ll for 228 galaxies in the UDF 
is available. While the spectroscopic sample does include 
28 redshifts from these data in the spectroscopic sample 
when based o n two emission lines with good signal to 
noise dMorris et al.ll2015l ) , the other 200 galaxy redshifts 
are not included because the 3D-HST redshifts are de¬ 
termined from the combination of the grism data and 
the photometric data. In essence this is a combination 
of grism emission or absorption lines and slope combined 
with the simultaneous determination of a photometric 
redshift with EAZY. While a pure photometric redshift 
from EAZY does not always agree with these redshifts, 
they are still based on the photometry and in such are not 
entirely independent tests for the photometric redshifts. 
In addition, while spectroscopic redshifts typically have 
a metric to select the more robust redshifts, the grism 
redshifts are presented without a quality flag. There¬ 
fore, we do not combine these redshifts with the robust 
spectroscopic sample defined in Section |4~T1 

While the 3D-HST grism redshifts are not independent 
of the photometry, they do provide another avenue to 
test the photometric redshifts. Since there are no qual¬ 
ity flags, the redshifts were manually inspected, and 214 
of the 228 redshifts were selected as robust (G. Bram- 
mer, private communication). Of these, 143 do not have 
robust spectroscopic redshifts and are a new test of the 
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Grism +EAZY redshift 

Figure 7. Redshift distribution of 214 grism redshifts. The or¬ 
ange represents all grism redshifts not covered by the spectroscopic 
redshifts, and the purple represents the 68 grism redshifts for which 
spectroscopic redshifts are available. 

photometric redshifts, and the distribution of the grism 
redshifts are shown in Figure [7] This includes 40 red¬ 
shifts at 2 > 2, almost twice the sample size of the 23 
spectroscopic redshifts at 2 > 2. 

As a test of the grism redshifts, the sub-sample of red¬ 
shifts that are also in our sample of robust spectroscopic 
redshifts are compared in Figure [5] This shows that of 
the 71 galaxies with both grism and spectroscopic red¬ 
shifts, three do not agree, and two of these would be an 
outlier by the definition used for the photometric red¬ 
shift comparison, | Az gr i sm \ /(I + z spec ) > 0.15. These 
two outlier redshifts are strongly influenced by a single 
emission line, while the other show no emissions lines at 
all, but the redshift is based on the SED and weak ab¬ 
sorption lines. While it is possible that the spectroscopic 
redshifts are in error, other work also find that the grism 
redshi fts someti me s disa gree with the spectroscopic red¬ 
shifts (jKriek et al.ll2014h . Additionally, the photometric 
redshifts agree with the spectroscopic redshifts, suggest¬ 
ing it is more likely that the grism redshifts are in error. 
Regardless, the agreement in Figure0is sufficiently good 
that it makes sense to compare the grism redshifts with 
the photometric redshifts, assuming that the grism red¬ 
shifts are correct. For the rest of the comparison with 
the grism redshifts, the three grism redshifts not in agree¬ 
ment with the spectroscopic redshifts are dropped from 
the analysis (e.g., Figure 0). 

The grism redshifts are compared to the photometric 
redshifts in Figure 0 The scatter is similar to that when 
comparing with the spectroscopic redshifts, ctnmad = 
0.031 and ctnmad = 0.040 for BPZ and EAZY. The OLFs 
are 2.9% and 4.7% for BPZ and EAZY, which become 
2.5% and 4.9% respectively for good ODDS and y 2 . The 
(Tnmad and OLFs remain the same when considering the 
sample with NUV data, as most of the grism redshifts fall 
in the area covered by NUV data. 

Similar to the spectroscopic sample, this is an improve¬ 
ment over the previou s cata logs when compared to the 
grism redshifts. The I Coe et all ; 200 lii redshifts have 
(tnmad = 0.031 and OLF o f 15% (12% with good ODDS 
and y 2 ), and the catalog bv lCameron et al.l (|2011h have 
(Tnmad = 0.059 and an OLF of 6.8%. This confirms the 
finding from the spectroscopic sample that the photo¬ 
metric redshifts presented here are robust, with signif¬ 
icant improvements over previous catalogs. There are 
also clear improvements when including the NUV data, 
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Figure 8. Grism versus photometric redshifts for the 71 galax¬ 
ies in the UDF with both grism and spectroscopic redshifts on a 
‘ pseudolog ’ scale (solid purple circles). An additional 6 galaxies are 
shown as open purple circles, whic h are from the sample of 14 grism 
redshifts excluded in Section 14.31 as potentially not being robust, 
and also have a spectroscopic redshift. The grism and spectroscopic 
redshifts are in good agreement except for three galaxies, showing 
that the grism redshifts are pretty good, although not as goo d th e 
sample of robust spectroscopic redshifts described in Section 14.11 
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Figure 9. Photometric versus grism redshifts for 211 galaxies in 
the UDF on a ‘pseudolog’ scale, showing the high quality of the 
photometric redshifts (excluding the three disagreements from the 
left panel plots). The blue circles represent redshifts for galaxies 
that are covered by the NUV data, and the green squares are for 
redshifts without NUV data. All redshift determinations include 
NIR and optical data. Filled symbols represent photometric red¬ 
shifts with good ODDS and x 2 ; ODDS > 0.9 and x 2 0( j < 4. The 
open symbols on the other hand are for redshifts That may have 
multiple peaks in P(z ) or may have a low quality fit, ODDS < 0.9 
or Xmod - > 4. Dashed lines represent outlier cutoff, defined as 
\Az\ /(I + Zspec) <0.15. 


especially in the OLFs as discussed below. 
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Figure 10. Photometric versus spectroscopic redshifts for galaxies in the UDF on a ‘pseudolog’ scale for the 125 galaxies with both 
spectroscopic coverage and NUV coverage. This shows the improved photometric redshifts with NUV and NIR data. Each panel is for 
redshifts determined by including different photometric bandpasses. The brown diamonds are redshifts with only the optical data used, 
the yellow triangles are redshifts including the optical and the NUV, and the green squares are redshifts including the optical and the NIR. 
Filled symbols represent photometric redshifts with good ODDS and Xmod > 4, and open symbols are for redshifts that may have multiple 
peaks in P(z) or may have a low quality fit, ODDS < 0.9 or Xmod > 4- Dashed lines in the bottom right panel is the outlier cutoff, defined 
as \A.z\ /(I + z S pec) < 0.15. 
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Figure 11. This figure is the same as Figure lTOl except including the 206 grism redshifts with NUV coverage instead of the spectroscopic 
sample. 
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4.4. Improvement of Photometric Redshifts with NJJV 

The improvements in the photometric redshifts with 
the addition of the NUV data can be considered by ex¬ 
amining the spectroscopic and grisrn subsamples that are 
covered by the NUV data (125 spectroscopic redshifts, 
206 grism redshifts). Before investigating this improve¬ 
ment, we note that the photometric redshifts presented 
here are already much improved over previous work due 
to the addition of the NIR data, even without the NUV 
data (similar to the redshifts determined_in the deep IR 
section of the UDF presented in iCameron et all (1201 ill , 
see Section 14.2.11) . However, the improvement of the 
redshifts by adding the NUV data is quantified by run¬ 
ning BPZ on the photometry without including the NUV 
data, and comparing to the spectroscopic and grism red¬ 
shifts. Removal of the NUV photometry results in an in¬ 
crease in ctnmad and OLF of the photometric redshifts. 
As mentioned before, for brevity only the single redshift 
code BPZ is considered. 

4.4.1. Spectroscopic and Grism Investigations 

For the spectroscopic redshift sample, the scatter im¬ 
proves slightly, from ctnmad = 0.029 without NUV pho¬ 
tometry, to onmad = 0.026 with NUV photometry. 
Moreover, the OLF decreases from 6.4% to 2.4% (no 
change with good ODDS and y 2 ), which is more than a 
factor of 2 decrease in the OLF when including NUV 
data. Note that the outlier fractions are susceptible to 
small number statistics, with the 6.4% corresponding to 
only 8 outliers. The photometric and spectroscopic red¬ 
shift comparison without NUV data are shown as green 
squares in Figure fTTTl 

For the grism redshift sample, the scatter again only 
improves mildly, from onmad = 0.034 without NUV 
photometry, to onmad = 0.031 with NUV photometry. 
However, the OLF decreases from 4.4% to 2.9% (4.0% 
to 2.5% with good ODDS and y 2 ), which again is close to 
a factor of 2 decrease in the OLF. The photometric and 
grism redshift comparison without NUV data are shown 
as green squares in Figure fill 

The improvement in the OLF with NUV data may be 
even more substantial than a factor of two, since the 
spectroscopic and grism redshifts preferentially sample 
z <2 (Figures |4] and [TJ). The NUV data are exp ected to 
aid considerably at z < 0.5 and 2 < z < 4 (e.g., iBenitezl 
120001: ITibert et al. 1120061 [Rafelski et aLll2009l) . due to color 
redshift degeneracies resolved by the NUV data. On the 
other hand, the NUV data do not significantly help the 
scatter in redshift, likely due to the other 8 photometric 
bands already sufficiently constraining the redshift, other 
than the possibility of outliers. 

In order to quantify the improvement in the redshifts 
with the addition of the NUV data further, the redshift 
determinations are compared for different combinations 
of bandpasses in Figures [TO] and eh The Figures show 
that using just the 4 optical bandpasses results in signifi¬ 
cant scatter and OLF, which is reduced by adding either 
NIR or NUV bandpasses. Quantitatively, for the spectro¬ 
scopic sample the 4 optical bands alone yield ctnmad = 
0.072, OLF=19.3%, the optical with 3 NUV bandpasses 
yield ctnmad = 0.071, OLF=13.0%, and the optical with 
the 4 NIR bandpasses yield ctnmad = 0.029, OLF=6.5%. 
This is compared to ctnmad = 0.026, OLF=2.4% for the 


spectroscopic sample including all 11 bandpasses. All 
outlier percentages are for good ODDS and y 2 on the sam¬ 
ple that is covered by NUV data. 

Similar statistics are obtained when examining the 
grism redshift sample, where the optical only bandpasses 
yield ctnmad = 0.064, OLF=18.4%, the optical and NUV 
bandpasses yield ctnmad = 0.060, OLF=14.0%, and 
the optical with NIR bandpasses yield ctnmad = 0.034, 
OLF=4.0%. This is compared to ctnmad = 0.031, 
OLF=2.5% for the grism sample including all 11 band- 
passes. 

This shows that either adding the NUV or the NIR 
significantly improves the redshifts, and including both 
yields the best results. The addition of the NIR data 
appears to provide a more significant improvement than 
the NUV overall, but not at all redshifts. For instance, 
the inclusion of the NUV data significantly reduces the 
scatter and OLF at z < 0.5 compared to the NIR data. 
This is accomplished by sampling the 4000A break with 
multiple filters, compared to only observing a flux decre¬ 
ment in the F435W band. 

Measurements of the Lyman break at z > 2 should 
also improve the redshifts of those galaxies, although the 
spectroscopic sample does not include sufficient galaxies 
at these redshifts to show this. While the grism sample 
includes a larger sample at z > 2, some fraction of the 
outliers observed in Figure EH are possibily due to the 
grism redshift being incorrect, similar to the three galax¬ 
ies in Figure [8] Also, the good agreement of the grism 
+ EAZY redshifts when including only the optical and 
NIR data are somewhat biased by the fact that the same 
photometry is already included in the determination of 
the comparison grism + EAZY redshifts. 

In addition, it is important to consider that the UDF 
NUV and NIR observations are not observed in an equal 
fashion. First, the NIR data have 4 bandpasses instead 
of the 3 for the NUV data. Second, the NIR data are 
deeper than the NUV data, consisting of 253 orbits of 
HST time compared to 46 for the NUV data. Third, the 
NIR filters are wider, have higher throughput than the 
NUV filters, and cover a larger wavelength range (Figure 
Lastly, most galaxies are significantly brighter in the 
observed NIR than observed NUV. All together, even if 
the NUV data improved the photometric redshifts more 
than the NIR, our data would not necessarily show it. 
Regardless, even in the presence of deep NIR data, the 
NUV reduce the OLF of the sample by a factor of 2. 

The scatter and OLF for any specific redshift range 
may be significantly more improved than for the entire 
sample. The NUV data is expected to im prove the color 
degeneracy at low a nd high redshifts (jllbert et al.1 120061 : 
iRafelski et al.ll2009th and indeed this is observed at low 
redshift, and insufficient spectroscopic redshifts exist at 
high redshift. The bottom left panel in Figure EOl which 
plots redshifts determined from optical and NIR data 
shows quite a bit of scatter at z < 0.5, while the top 
right panel for the redshifts determined with optical and 
NUV data shows very little scatter at those redshifts. 

4.4.2. Comparison of Redshifts with and without NUV 

The photometric redshifts obtained with and without 
NUV data are compared to each other in Figure EH to 
show how the NUV data alters the redshifts for the 6459 
galaxies brighter than F606W magnitude of 29 and with 
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Xmod < 4. This magnitude cut is based on the depth 
of the NUV data and the typical colors of galaxies in 
the sample. Galaxies brighter than this magnitude may 
have their redshifts improved with the addition of the 
NUV data, while fainter galaxies are unlikely to show 
improvements. Together, the magnitude and X^ od cuts 
result in a sample of good redshift fits in which the NUV 
data may contribute. 

There are two scenarios that are evident in Figure fl2l 
In the first, the photometric redshifts without the NUV 
data have 0DDS> 0.9, and therefore a change in the pho¬ 
tometric redshift with the addition of NUV data is not 
expected, as the probability function suggests only one 
possible redshift for that photometric fit. The second is if 
the redshifts have 0DDS< 0.9, in which case the redshifts 
are uncertain and a change in redshift is not unexpected. 

An outlier fraction can be defined again by assuming 
that the photometric redshift with NUV data is correct, 
and requiring \&.z p hot\ /(I + z spec ) < 0.15. In this case, 
7.0% of the eight band optical and infrared redshifts 
would be at the wrong redshift. On the other hand, 
only 2.4% do so with good ODDS. If the redshifts deter¬ 
mined with the NUV data are correct, then the inclusion 
of the NUV data removes at least 2.4% of the outliers 
with good ODDS and Xmod- This 2S a similar to the 
improvement observed in the OLF for the spectroscopic 
and grism samples when adding NUV data, suggesting 
that the outlier rate is consistent with an improvement 
by a factor of 2 when including NUV data as suggested 
in the much smaller spectroscopic and grism samples. 

4.4.3. Number of Galaxies with Single Peak in P(z) 

Another method to estimate the improvement of the 
photometric redshifts with the addition of the NUV data 
is to consider the increase in the number of galaxies with 
a single peak in P(z) (which therefore have unambiguous 
redshifts) with the addition of the NUV data, as done in 
the predictions in iTeplitz et all (12013( 1 . This is similar 
to requiring a good ODDS value, but provides a little 
more information on P(z ), as it is possible for a redshift 
to have 0DDS> 0.9 and still have a second peak in the 
probability distribution. 



Figure 12. Comparison of photometric redshifts with and with¬ 
out including the NUV data in the photometric fits for the 6459 
galaxies with F606W magnitude < 29 and with Xmod < T The 
pink circles have 0DDS> 0.9, and the gray squares have 0DDS< 0.9. 
This shows that the addition of NUV data significantly affects the 
redshifts of galaxies in the UDF. If the redshifts with NUV data 
are correct, then the NUV data removes 2.2% of the outliers with 
good ODDS and x^ od - 


Figure [13] shows the number of galaxies with a single 
peak as a function of redshift for different numbers of 
bandpasses included in the photometric fit. The photo¬ 
metric redshifts improve with the addition of each set 
of bandpasses, with the inclusion of all the bandpasses 
yielding the largest number of unambiguous galaxy red¬ 
shifts. In the redshift range 1 1$ z 3 either adding the 
NUV or the NIR data yield similar improvements, while 
the addition of the NIR improves the redshifts more at 
1 1$ z 1$ 2 and the NUV improves the redshifts more at 
z > 2. This is sensible, as the addition of the NIR data 
samples the 4000A break at lower redshifts, and the addi¬ 
tion of the NUV data enables sampling the Lyman break 
for galaxies at higher redshifts. While including both the 
NIR and the NUV data increases the number of single 
peaked galaxies the most, the improvement from adding 
either to the other two is not as significant as with either 
the addition of the NIR or NUV to just the optical. 

The increase in the number of galaxies with a single 
peak in P(z) with the addition of the NUV data to the 
optic al and NIR data is less than the number presented in 
(ITeplitz et alJl2013h . although the total number of galax¬ 
ies is similar. There are multiple factors contributing to 
this difference. First, the number of galaxies predicted to 
have good red shifts with op t ical a nd NIR data is larger 
than shown in ITeplitz et all (1201311 , which results in less 
change with the addition of t he NUV da ta. The cause of 
the lower numbers in ITeplitz et al.l ( 2 0131 1 is likely related 
to an increased number of filters and sensitivities from 
the UDF 12 data that were not included in this o riginal 
estimate dKoekemoer et al.ll2013HEllis et aJJl2013f) . 

Second, the NUV data presented h ere are less se nsitive 
than those used in the predictions in lTeplitz et al.1 (l2013l l 
because this paper only uses the half of the data which 
are not binned and include the addition of postflash to 
the images (making them more robust to CTE degra¬ 
dation, although also less sensitive). The other more 
sensitive half of the data is binned, and no pixel-based 
CTE corrections exist yet for binned data. Since the pho¬ 
tometry without such corrections would be incorrect (see 
ITeplitz et al.ll2013l) . those data are not included here. If a 
pixel-based CTE correction for the binned data could be 
made, it would improve the redshifts further, especially 
for fainter galaxies. 

5. CATALOG OF THE UDF 

The UDF catalog is made available in both a FITS ta¬ 
ble and an ASCII table at the UVUDF website^ an d on 
MASlFl As the table has a large number of columns, the 
FITS table is recommended for use. The catalog columns 
are described in detail in Table [5] The final catalog is 
trimmed to only include sources that fall within the cen¬ 
tral 11.4 arcmin 2 of the image based on a minimum 30% 
exposure time in the optical UDF. This cut reduces the 
area covered by the shallow NIR data to 6.8 arcmin 2 . 
However, measurements are made in the full UDF field 
of view for all bandpasses, enabling reliable photometry 
of sources near this artificial edge, but objects on the ac¬ 
tual UDF edge are not included. In addition, the catalog 
is trimmed to include only sources with two photometric 

20 http://www.rafelski.com/uvudf/catalogs.html 

21 http://archive.stsci.edu/prepds/uvudf/ 
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Figure 13. Number of galaxies with unambiguous photometric redshifts (a single peak in their P(z)) as a function of redshift for 
photometric redshifts including different numbers of bandpasses. This shows how different bandpasses help improve the photometric 
redshifts over the base redshifts determined from only optical data (brown). Photometric redshifts including only the NUV and optical 
bandpasses are in gold, those including only the NIR and optical are in green, and those including the NUV, optical, and NIR bandpasses 
are in blue. The left panel shows the improvement obtained by including the NUV and optical data (gold), and the right panel shows the 
same for the NIR and optical data in green. The NIR, improves the redshifts more at 1 < z < 2 and the NUV improves the redshifts more 
at z > 2, as the bandpasses sample the 4000A and Lyman break respectively. 


measurements at > 3tr confidence. While a small num¬ 
ber spurious sources could still be present in the catalog, 
this removes the majority of them. This also removes 
high redshift sources that would only be detected in the 
F160W bandpass. 

The final catalog ID numbers identify which of the 
four SExtractor runs each source measurement is from, 
as described in Section 13.51 Sources from the deep- 
detection and normal deblending thresholds have ID 
numbers in the range 1-19,000, while those with shal¬ 
low detection and normal deblending range from 20,GOO- 
29,000. Sources from the normal detection and low de¬ 
blending have ID numbers in the range 30,000-49,000, 
while those from the shallow detection and low deblend¬ 
ing range fro m 50 , 000- 59,000. The corresponding ID 
numbers from ICoe et all (120061) are included for back¬ 
wards compatibility, matched using the seg mentation 
maps of both this catalog and that of lCoe et al.l (|2006f ). 
However, note that any single source m ay not be matched 
to an object in the ICoe et ahl (120061) catalog, and sirni- 
larly , two differe nt sources may be matched to the same 
ICoe et all (120061 ) source. 

The photometric magnitude and flux values are pro¬ 
vided for each filter. The total magnitudes are aperture- 
matched PSF corrected total measurements, including 
corrections for Galactic extinction. These total measure¬ 
ments are based on the isophotal color and the total flux 
measured via the Kron aperture flux (see Section EU). 
The aperture correction to convert from the isophotal 
magnitude to total magnitude is available in the cata¬ 
log, along with the original isophotal flux values. The 
NUV photometric measurements and resulting photo¬ 
metric redshifts all are based on the smaller F435W aper¬ 
tures as described in Section lTBl The table also indicates 
whether a source is covered by the NUV data, and if so, 
if it falls near the NUV edge or chip gap. The optical 
and NIR data cover the entire UDF FOV, and therefore 
no flag is provided. In addition, some basic morphologi¬ 
cal measurements from the F775W SExtractor runs are 
also provided, such as FWHM and Ellipticity. 

The photometric redshifts from both BPZ and EAZY 
are presented, along with the number of filters available 
to the codes, and the number of measurements above 5a 


included in the fits. The quality of the photometric red¬ 
shifts are included, described by ODDS, Xmo<u an d xt as 
described in Section [4j Care should be taken when us¬ 
ing these redshifts to select galaxies with the appropriate 
values; to obtain a robust sample, we recommend using 
a minimum ODDS of 0.9, maximum Xmod °f 4 for BPZ 
redshifts, and a maximum xt °f 10 for EAZY redshifts. 
Stricter cuts may be used for a higher confidence pho¬ 
tometric redshift, at the cost of a smaller sample size. 
Spectroscopic and Grism redshifts are also provided if 
available, and known stars are marked. 


6. SUMMARY 


New NUV and NIR imaging obtained with WFC3 on 
HST provides the opportunity to significantly improve 
the photometric redshifts of galaxies in the UDF. The 
NUV data come from the UVUDF imaging campaign 
including three filt ers, F225W, F275W, and F336W 
(iTeolitz et al.l [2013 ). and cover ~ 60% of the UDF 
FOV. The NIR data is obtained from a combination of 
the UDF09, UDF12, and CANDELs data in order to 
obtain the deepest NIR coverage and c over the entire 
FOV of the UDF (lOesch et al.l 2010dlbl: iBouwens et ahl 


2011bl: iGrogin et alJl2011t iKoekemoer et all 1201 ll 120131 : 
Ellis et. a, 1.112013D . 


The NUV data are newly calibrated as described in 
Section [2] and the Appendix, and made available to the 
public through MAST. In short, the data are corrected 
for CTE degradation, and calibrated with custom dark 
Hies, as the current STScI CDBS (CRDS) dark calibra¬ 
tion files are no longer sufficient to calibrate the NUV 
data largely due to CTE degradation of the detector. 
The new custom dark files are CTE corrected, mostly 
remove a gradient and blotchy pattern from the science 
data, and properly flag hot pixels. The custom process¬ 
ing significantly improves the final mosaic image quality. 
Darks processed in a similar fashion should be used for 
any observations with low backgrounds, and STSCI plans 
to release such darks sometime in 2015. In order to min¬ 
imize the blotchy pattern in future data, we recommend 
that observations should use dither patterns larger than 
~ 20 pixels to reduce the pattern in the final mosaics. 

Sources are detected in an eight-band-averaged image 
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Table 5 

Column Description of UDF Catalog 


Column No. 

Column Title 

Description 

i 

ID 

Object identification number 

2 

COEJD 

Obiect identification number from Coe et al. (2006) a (-99 if nojriatch) 

3 

RA 

Right ascension (J2000 in units of decimal degrees) 

4 

DEC 

Declination (J2000 in units of decimal degrees) 

5 

X 

Image X pixel coordinate in the UDF mosaic 

6 

Y 

Image Y pixel coordinate in the UDF mosaic 

7-17 

MAG.* 

Total AB magnitude of each filter b 

18-28 

MAGERR.* 

Total AB magnitude uncertainty of each filter 0 

29-39 

FLUX.* 

Total flux of each filter in units of ^Jy d 

40-50 

FLUXERR.* 

Total flux uncertainty of each filter in units of ^Jy d 

51-61 

FLUX.ISO.* 

SExtractor isophotal Flux of each filter in units of e — /s d 

62-72 

FLUXERR.ISO.* 

SExtractor isophotal Flux uncertainty of each filter in units of e — /s d 

73 

APCOR 

Aperture correction to convert from isophotal magnitude to total magnitude 6 

74 

FWHM 

SExtractor F775W FWHM in units of pixel (1 pixel = 0.03 // ) 

75 

AREAF 

SExtractor F775W isophotal area (filtered) above detection threshold in units of pixel 2 

76 

STELLARITY 

SExtractor F775W stellarity 

77 

ELLIPTICITY 

SExtractor F775W ellipticity 

78 

THETA 

SExtractor position angle in units of degrees 

79 

UVUDF.COVERAGE 

1 = covered by NUV data, 0 = not covered by NUV data 

80 

UVUDFJ3DGEFLG 

1 = close to NUV edge or chip gap, 0 = not close to edge or chip gap 

81 

Z.BPZ 

Bayesian photometric redshift (BPZ) 

82 

ZMIN.BPZ 

BPZ 95% lower limit to redshift 

83 

ZMAX.BPZ 

BPZ 95% upper limit to redshift 

84 

ODDS.BPZ 

BPZ integrated P(z) contained within 2 * 0.03 * (1 + z) 

85 

CHISQ2.BPZ 

BPZ modified reduced chi square (Xmod) 

86 

TEMPLATE.BPZ 

BPZ template number f 

87 

Z.EAZY 

EAZY photometric redshift 

88 

ZMIN.EAZY 

EAZY 95% lower limit to redshift 

89 

ZMAX.EAZY 

EAZY 95% upper limit to redshift 

90 

ODDS.EAZY 

EAZY integrated P(z) contained within 0.2(1 + z) 

91 

CHISQ.EAZY 

EAZY reduced chi square (x2) 

92 

NFOBS 

Number of filters available for photometric redshift 

93 

NF5SIG 

Number of filters with signal to noise above 5 

94 

SPECZ 

Spectroscopic redshift (-99 if no value) 

95 

SPECZ.REF 

Reference for spectroscopic redshift 8 

96 

GRISMZ 

Grism + EAZY redshift from 3D-HST team (Brammer et al. 2012) (-99 if no value) 

97 

STAR 

Stars identified by GRAPES program (Pirzkal et al. 2005) (l=star) 


Note. — Column information for the UDF catalog, available both as a fits and ASCII file. Multiple entries are for all 11 bands, where the * 
represents each filter described in Table |T] F225W, F275W, F366W, F435W, F606W, F775W, F850LP, F105W, F125W, F140W, and F160W. 
a Matched based on segmentation map overlap. (-99 if no match) 

k The total AB magnitude based on the Kron radius, including extinction, aperture, and PSF corrections. 99 if not detected at la. -99 if not 
covered by a filter, or for NUV filters, also -99 if not detected in F435W image. 

c If not detected then equal to the lcr limiting AB magnitude. -99 if not covered by a filter, or for NUV filters, also -99 if not detected in 
F435W image. 

d If not detected then still provides measured flux or flux uncertainty. -99 if not covered by a filter, or for NUV filters, also -99 if not detected 
in F435W image. 

e Includes extinction and PSF correction 

^ BPZ templates for 11 galaxies as shown in Figure[3] (1) E117 (2) E116 (3) E115 (4) E114 (5) ESO (6) Sbc (7) Scd (8) SB1 (9) SB2 (10)SB3 
(11) SB11. The 9 interpolated galaxies between adjacent galaxies are represented as decimal places between the integer galaxy numbers. 

8 (l)[Lc Fcvrc ct al. (2004) (2) Szokolv ct al. (2004) (3) Mignoli ct al. (2005) (4'Daddi ct al. ( 20051) (5) IVanzella et all (|20 05. 2006, 2008, 2009) 
(6) IPopesso et al ] l12H09li (7) IBalestr a ct al. (2010) (8) Kurk ct ^ (1201311 (9) IMorris et all d2015fl 


including the optical and NIR. Aperture-matched PSF 
corrected photometry is performed on the eleven photo¬ 
metric bandpasses, and are included in Table [5] Numer¬ 
ous SExtractor runs per filter were used to optimize the 
detection and segmentation definitions of sources in the 
UDF. Specifically, two different detection thresholds and 
two different deblending thresholds were run, and then 
merged together. In addition, the NUV data use sep¬ 
arate aperture definitions determined from the F435W 
image to maximize SN, and minimize systematics from 
any leftover blotchy pattern in the NUV data (See sec¬ 
tion [A3]). 

Photometric redshifts were derived using two dif¬ 
ferent redshift codes, B PZ (iBenfted I2000H and 
EAZYlBr ammer et all 13)081 ) using different galaxy tem¬ 


plates and priors. These redshifts are compared to a 
new compilation of 169 reliable spectroscopic redshifts of 
galaxies (Table [4|. This comparison reveals a low scat¬ 
ter of unmad = 0-028 for BPZ and ctnmad = 0.030 for 
EAZY, and low outlier fractions (OLF) of 2.4% and 3.8% 
for BPZ and EAZY respectively for good ODDS and y 2 . 
Results from the comparison to grism redshifts are sim¬ 
ilar, with (Tnmad = 0.031 and onmad = 0.040 for BPZ 
and EAZY and OLF of 2.5% and 4.9% respectively for 
good ODDS and y 2 . This is an improvement of ~ 2 in 
ctnmad and a factor of ~ 3 in OLF o ver t he l ast co mpre- 
hensive UDF redshift catalog bv lCoe et al.i (l2006i l. 

We showed that adding the NUV data to the photo¬ 
metric redshift derivations in addition to the optical and 
NIR gave a mild improvement in ctnmad, an d a factor 
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of ~ 2 improvement in the OLF. The improvement of 
the redshifts with adding NUV or NIR data to the opti¬ 
cal data depend on the redshift, with the NIR improving 
the redshifts at 1 < z < 2 more, and the NUV doing so 
at z < 0.5 and z> 2. In addition, the NUV data appear 
to significantly reduce the scatter and OLF at 2 < 0.5. 
It is important to consider the caveat that in these com¬ 
parisons, the NUV data only consist of 46 orbits over 3 
filters, versus the 253 orbits of NIR data over 4 filters. In 
other words, for a smaller investment of NUV observing 
time, the outlier fraction can be significantly improved 
with the addition of the NUV. Such improvements in the 
photometric redshifts are observed even when consider¬ 
ing only the NUV and optical data, showing the power 
of including the NUV data. 

The photometry and photometric redshifts of all the 
sources in the full UDF are presented in Table [5] as both 
a FITS and ASCII table. Overall, this catalog provides 
photometry in a new wavelength regime, and signifi¬ 
cantly improves the photometric redshifts. It will aide 
future galaxy research in the UDF, and has already con¬ 
tributed to a n umber of studies using the NUV data and 
redshifts (e.».. [Rond et al.ll20l4 iKurczvnski et al.ll20l4 
iMei et alJl2014l ) and others in preparation. 
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A. WFC3/UVIS DARK CALIBRATIONS 

The UVUDF data analysis revealed a number of im¬ 
provements that can be made in the creation of the dark 
calibration files. We worked closely with STScI to solve 
problems with the darks and communicate the best prac¬ 
tices for future observing strategy and data reduction. 
The need for improved darks are caused by changes in the 
characteristics of the detector (such as CTE degradation) 
since the methodology for dark creation for WFC3/UVIS 
was developed. There are three effects discussed here: 
hot pixel detection, background gradients, and a blotchy 
background pattern. The official dark calibrations are 
part of the Calibration Reference Data System (CDBS), 
and the STScI released darks will hereafter be referred to 
as CDBS darks. We note that recently STScI switched 
over to the Calibration Reference Data System (CRDS), 
which currently uses the same CDBS darks. 

Part of the WFC3/UVIS calibration process is the sub¬ 
traction of a dark reference file to remove the dark cur¬ 
rent and to identify hot pixels that can cause significant 
artifacts in the images. STScI releases new darks every 


3-5 days, which is necessary due to the appearance of 
~ 500 new hot pixels per day. The increase in hot pixels 
is mitigated by annealing the detect or once per month, 
removing > 70% of the hot pixels ({Borders fc Baggettl 
12009 1. Even so, the number of permanent hot pixels 
which are not removed by anneals is increasing by 0.05- 
1% per month (WFC3 instrument handbook). 

The current methodology used for the CDBS darks is 
to combine the raw darks obtained over a 3-5 day period 
to create superdarks for each time period, composed of 
~ 10—20 dark exposures with integration ti mes of ~ 900s 
each. The process is described i n deta il in iMartel et al.1 
(|2008f) and IBorders fe Baggettl (|2009f ). but the general 
process is outlined here. First, the process cleans the 
raw darks of cosmic rays, and then creates an average of 
all the raw dark files. A threshold is then used to find 
all pixels deviant from the median dark value. These 
pixels are marked as hot, and then the final superdark 
is the median value of the average dark (a single value 
for all good pixels), with the hot pixels superimposed 
and marked in the data quality array. This is done to 
minimize introducing noise in the data from the uncer¬ 
tainty in the dark current per pixel from a small number 
of exposures over the 3-5 day period. While this process 
was sufficient for the data products shortly after instal¬ 
lation on HST, the detector characteristics have changed 
making this process insufficient for current WFC3/UVIS 
images. 

A.l. Missed Hot Pixels 

The first issue with the CDBS darks is that an out¬ 
dated def inition of a hot pixel is ap plied (originally deter¬ 
mined in iBorders fe Baggettl 1200911 . which results in un¬ 
masked warm-to-hot pixels remaining in data obtained 
in recent years. Under this definition, pixels with val¬ 
ues > 0.015e _ /s are flagged as hot. However, since the 
time that this threshold was determined the characteris¬ 
tics of the detector have changed, especially due to CTE 
degradation. The situation is worse because the CDBS 
darks are not CTE corrected, resulting in hot pixels being 
missed even if they would otherwise have been above the 
threshold level, and are present in the CTE corrected sci¬ 
ence data. The leftover hot pixels that remain unmasked 
yield significant artifacts for mosaics calibrated with the 
CDBS darks. 

These missed hot pixels are more prevalent far from the 
readout, where CTE degradation plays a more important 
role. This effect is evident in Figure Qj] which shows the 
number of hot pixels as a function of row number, where 
the center of the two chips is at pixel row 0, and the 
readout of the two CCDs are at the left and right sides 
of the plot. The black line shows the number of hot 
pixels per row in the first in-flight CDBS dark released 
in 2009, and a constant number of hot pixels per row is 
observed. All lines on this plot after this are expected to 
have an increased number of hot pixels with time, due 
to the growth of the nu mber of permanent hot pixels 
(IBorders fc Baggettll200 9). However, the gray line shows 
the hot pixels for one of the dark files associated with 
the UVUDF data in 2012. At the edges of the plot, the 
number of hot pixels is indeed higher than the black line, 
but at the center of the plot far from the readouts, the 
number of hot pixels per row decreases significantly. 

There is no reason to expect that the number of hot 
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Pixel Row Number 

Figure 14. The number of hot pixels per row versus row number 
for WFC3/UVIS data determined from dark data. The center of 
the two CCDs is at pixel row number 0 at the center of the plot, 
and the edges of the chips near the readouts are at the right and left 
side of the figure. The black, gray, and yellow lines are determined 
from CDBS darks from 2009, 2012, and 2014 respectively, with the 
2014 dark based on raw darks that are post-flashed. The green 
and blue lines include the more aggressive hot p ixel m asking from 
the new dark methodology described in Section roi for the same 
2012 and 2014 darks as the previous ones. The pink, green, blue, 
and brown lines also include an additional CTE correction step on 
the raw dark files, before the creation of a superdark. The brown 
line is for a 2014 superdark obtained with high level of postflash to 
minimize CTE trails. As expected, the decrement of the number 
of hot pixels in rows near the center of the chip is not present in 
this dataset. This figure shows that the CDBS darks from 2012 
and 2014 are missing a substantial number of hot pixels, especially 
for rows far from the readout. The yellow and blue lines also shows 
the improvement in hot pixel detection in the post-flash data. 

pixels is physically developing differently on different 
parts of the chips. This effect is due to the CTE degra¬ 
dation causing hot pixels to be missed far from the read¬ 
out. Therefore, a significant number of hot pixels are 
not masked in the CDBS darks at later times, as the 
gray line would be expected to be flat from one side of 
Figure Il4l to the other. Therefore, the 2012 CDBS darks 
are missing 5 7% of the ho_t_pixels deemed important to 
mask by ([Borders fc Bagget5l200 9b 

The missed hot pixel issue is somewhat reduced by the 
introduction of post-flash darks in October 2012, and 
the hot pixels per row in a dark from 2014 is shown in 
yellow in Figure [14] However, even in the post flash 
data, 35% of the hot pixels are missed. If these data 
are CTE corrected before being processed in otherwise 
the same manner as the CDBS darks (pink line), then 
25% of the hot pixels are missed. This shows that even 
if a CTE correction is applied to the post-flashed raw 
darks, there are still a lot of hot pixels missed. Since 
the CTE correction code aims to not over correct the 
pixels, it is likely that these hot pixels are missed due to 
imperfect CTE corrections. The fact that the post flash 
darks (yellow and pink curves) are flatter than the darks 
without post flash (gray curve) in Figure [TJ] is evidence 
that a large fraction of the missed pixels are due to CTE 
degradation. 

A.2. Background Gradient 

The second issue with the CDBS darks is that the me¬ 
dian value of the average darks is applied as the value 
of all pixels in the dark frame. This median dark file is 
not suitable for low background data, because it leaves a 
low-level gradient. This gradient is typically small com¬ 
pared to the sky background in exposures with high sky 



Pixel Row Number 

Figure 15. The average background flux per row versus row 
number for WFC3/UVIS darks. The center of the two CCDs is 
at pixel row number 0 at the center of the plot, and the edges 
of the chips near the readouts are at the right and left side of 
the figure. The black, red, green and blue lines are from darks 
obtained in 2009, 2012, and 2014. The green, blue, and orange 
lines also include a CTE correction. The orange line represents 
the dark background used for calibrating the post-flashed UVUDF 
data based on early test po stflas hed dark exposures obtained in 
2012 as described in Section The CTE correction and post¬ 
flashing the raw darks both significantly reduce the background 
gradient, but does not remove it. 

backgrounds, however, in the low-background UV imag¬ 
ing, it is the dominant structure. The use of CDBS darks 
therefore results in a significant gradient in the science 
mosaics. 

The background gradients measured in the averaged 
darks for different time periods are shown in Figure 1151 
Similar to Figure [ll] the center of the two CCDs is at 
pixel row number 0 at the center of the plot, and the 
edges of the chips near the readouts are at the right and 
left side of the figure. The black line shows the first in¬ 
flight dark background gradient, and it is relatively flat. 
This explains why the CDBS darks do not attempt to 
correct this background gradient. However, the darks 
associated with the UVUDF in 2012 have a significant 
gradient, shown as the pink line in Figure [15] The back¬ 
ground level is significantly higher at the center of the 
chips far from the readout than at the edges close to the 
readout, with a factor of four variation in the background 
level from the center to the edge of the chip. This sug¬ 
gests that CTE degradation may be responsible for the 
observed gradient, with flux from cosmic rays and hot 
pixels smeared out into the background. 

The green line in the plot shows the gradient from the 
same data as the pink line, but with a CTE correction 
applied. This dark shows a reduction in the background 
gradient, confirming that at least part of the gradient 
from the pink line is indeed due to CTE degradation. In 
addition, the blue line shows a dark obtained in 2014 with 
post-flash enabled and with a CTE correction applied. 
Although the gradient is further reduced, it still has a 
factor of two variation across the chip. 

The most likely cause is that the CTE degradation 
strategies and corrections do not fully correct the data 
for the CTE degradation, although there could be other 
physical causes for the remaining gradient in the back¬ 
ground. One such possibility is if the postflash light is 
not fully subtracted, as the postflash illumination is not 
flat, but i nstead is ~ 3 0% b righte r near the center of 
the chip (|MacKentv fc Smithll2012if . A slight imperfec- 
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Figure 16. Both panels show the same empty region in the UDF 
in the F275W filter with WFC3/UVIS. The top panel is processed 
with the CDBS darks released by STScI, and the bott om panel 
with the new improved darks described in Section IA.4I The top 
panel exhibits a blotchy pattern which is significantly reduced in 
the bottom panel after applying the improved darks to the science 
data. 


tion in the subtraction of the postflash could contribute 
to the ~ le _ gradient left over in the dark observations. 
Since the science data are flashed at the same level as the 
darks, any such light would be removed when the darks 
are subtracted and thus could be missed in tests. 

A.3. Blotchy Background Pattern 

In addition to the background gradient, the dark back¬ 
ground exhibits a blotchy pattern on scales of ~ 20 pix¬ 
els. This pattern is difficult to characterize, as the signal 
to noise (SN) in each 3-5 day dark is very low. In ad¬ 
dition, the darks contain many hot pixels and artifacts, 
which make the noise pattern difficult to visualize in the 
darks themselves. Therefore, the pattern is character¬ 
ized and shown in the science data in an empty region 
of the UDF. The top panel in Figure El shows an empty 
region of the UDF with a non-uniform background noise 
pattern. This background pattern exhibits blotches that 
are ~ 0.001 — 0.002 e _ /s in magnitude (both positive 
and negative) with ~ 20 pixel diameters (~ 0.8 arcsec- 
onds). This translates to changes of ~ 1 — 2 e _ in the 
original raw dark frames. Given this blotchy pattern in 
the data, future observations should use dither patterns 
larger than ~ 20 pixels to reduce the pattern in the final 
mosaics. 


A.4. New Dark Methodology 

Due to these limitations of the CDBS darks, and the 
need for cleaner image mosaics in the deep UDF obser¬ 
vations, we developed a new methodology to create im¬ 
proved darks for WFC3/UVIS superdarks that mitigate 
the issues described above. The new dark methodology 
begins by applying the pixel-based CTE correction on 
all the individual raw dark files. Those data are then 
used in two different (but connected) steps; the first is 
to determine the hot pixels, and the second is to pro¬ 
duce the best background level. For the first step, the 
darks from a 3-5 day window are used to create a su¬ 
perdark in a somewhat similar manner as done by STScI 
in creating the CDBS darks. The same time frame as 
the CDBS darks is employed, and thus the data are on 
the same cadence as CDBS. A short cadence is necessary 
for masking the hot pixels, since they vary rapidly. The 
hot pixels are determined in a different fashion than the 
CDBS darks. For each 3-5 day dark, the background per 
averaged dark is modeled with a 3rd order polynomial, 
which is used to remove the gradient. After the gradient 
is removed, hot pixels are found with a more aggressive 
threshold. This step is used to define the hot pixels in 
the darks on short time scales, and also provides cosmic 
ray cleaned raw darks for the later step. 

A.4.1. Finding Hot Pixels 

The threshold for finding hot pixels is set based on 
the darks in 2012 and has been optimized for the UDF 
observations. A different threshold is used for the post- 
flashed and non-postflashed darks. The thresholds are 
set such that the number of hot pixels per row at the 
center of the chip is equal to the number of hot pixels 
in the original CDBS darks close to the readout. In this 
fashion, all hot pixels that would have been originally 
masked are found. This results in a 3.7cr threshold for 
non-postflashed data, and a 4.9er threshold for post-flash 
data. The standard deviation of the pixels is determined 
from an iterative 3er rejection, thereby removing any hot 
pixels from the distribution before a is determined. Also, 
the background gradient is removed before applying the 
threshold, such that the gradient does not cause different 
numbers of hot pixels to be found across the chip. We 
ran tests with different thresholds to determine if a higher 
or lower threshold would produce better science images, 
and found that thresholds similar to those determined 
by STScI produced the right balance of masking all the 
noticeable hot pixels, while not masking too much of each 
raw data frame. 

The resulting number of hot pixels as a function of 
pixel row number for our new darks are shown in Figure 
El The green line shows the resultant hot pixel numbers 
for a CTE corrected dark from 2012, and the blue line 
shows the same for a 2014 CTE corrected post-flashed 
dark. The number of pixels at the center (pixel row 0) 
in the green and blue lines match the number of hot 
pixels per row near the readout for the equivalent CDBS 
dark as shown in the gray and yellow lines respectively, 
confirming the threshold levels. In order to mask the hot 
pixels far from the readout, a higher fraction of pixels 
are masked than actual hot pixels in the data. 

As a test of the CTE degradation affecting hot pix¬ 
els, a recent set of calibration darks obtained by STScI 
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in March 2014 were post-flashed to a very high level, 
~ 90e _ . This enables a check on the cause of the vari¬ 
ation in number of hot pixels per row across the chip, 
as darks with a very large postflash count are expected 
to experience fewer charge traps, and therefore less CTE 
loss. The brown line in Figure Q3] is from a superdark 
processed using the methodology described above, and 
is mostly flat, which is consistent with the expectation 
that when CTE losses are minimized, hot pixels are not 
preferentially lost far from the readout. However, the 
postflash subtraction with such a large postflash is not 
perfect, and the darks contain residual postflash light. 
Therefore, investigations of the background gradient is 
not possible with this dataset, nor can one apply the 
original CDBS threshold for a hot pixel. The slight in¬ 
crease in the number of hot pixels at the center of the chip 
is likely due to the ~ 30% variation in the illumination 
pattern of the postflash, changing the noise properties at 
the center of the chip where the illumination is brightest. 

One downside of more aggressive hot pixel flagging is 
that a lot of pixels are flagged as hot that are likely just 
warm pixels, and do not need to be flagged as hot. This is 
especially evident in Figure [14] for rows near the readout, 
where a factor of more than two times the number of hot 
pixels are flagged than the rows furthest away from the 
readout. However, the fraction of unnecessarily flagged 
pixels of the detector is small, at ~ 1% and ~ 2% of 
the chips for data with and without post-flash respec¬ 
tively. Even with this increase, the total fraction of hot 
pixels flagged remains small, at ~ 3% and ~ 4% for with 
and without post-flash respectively. On the other hand, 
if the same thresholds are applied to the initial inflight 
2009 darks (black line), the number of hot pixels would 
be doubled. Regardless, this increased hot pixel flag¬ 
ging should not be of major consequence, since the total 
number of hot pixels remains small, and flagging the hot 
pixels improves image quality sufficiently to justify this 
choice. 

A.4.2. Removing Gradient and Blotchy Pattern 

The second step in the new dark methodology is to ob¬ 
tain the best dark background to handle the gradient and 
blotchy pattern described above. Specifically, we create 
a dark with spatial structure in it, rather than a dark 
with a single value for all good pixels. To do so, the 
CTE corrected cosmic ray cleaned images from a single 
anneal cycle are averaged together, with the hot pixels 
from each 3-5 day window being masked. This deter¬ 
mines the actual dark level for each good pixel. This 
averaged dark per anneal cycle is then used in conjunc¬ 
tion with the hot pixel map on the 3-5 day period to 
create a new superdark consisting of the average dark 
over the entire anneal cycle, with hot pixels masked at 
the shorter cadence. We determined that using a shorter 
time period than an entire anneal cycle did not improve 
the blotchy pattern, and the characteristics of the darks 
changed sufficiently between anneals to preclude averag¬ 
ing together raw darks from different anneal cycles. Also, 
note that the modeled dark gradients in the first step are 
not used other than to find the hot pixels. 

The average dark per anneal intrinsically includes any 
gradient in the dark background, as well as the observed 
blotchy pattern. The resulting darks significantly im¬ 
prove the image quality in the UVUDF, as well as in 


other WFC3/UVIS programs. The improvement is evi¬ 
dent in Figure 1161 where the top panel shows an empty 
region of the NUV UDF data processed with the CDBS 
darks, and the bottom panel shows the same data pro¬ 
cessed with the new dark calibrations. While the blotchy 
pattern is not removed completely, it is significantly re¬ 
duced. 

In order to obtain the best hot pixel mask and dark 
background, it is best to match post-flashed science data 
with post-flashed darks, and CTE corrected science data 
with CTE corrected darks. Not doing so can yield missed 
hot pixels in the science frames, and imperfect dark back¬ 
ground levels. The UVUDF data were obtained before 
STScI started obtaining post-flashed darks on a regular 
basis. However, 30 dark frames with post-flash were ob¬ 
tained during the testing process in the same month as 
the UVUDF science images, and are the final dark expo¬ 
sures used to calibrate the UVUDF. 

A comparison of the science mosaics shows that the 
blotchy pattern is less prominent when the science data 
are processed with darks based on these 30 darks, rather 
than the higher SN darks produced from a larger num¬ 
ber of non-postflashed darks obtained at the same time. 
Since no post-flash darks exist on the same days as 
the UVUDF science exposures, we include the hot pixel 
masks from the non-postflashed darks in addition to the 
hot pixel mask from the post-flashed darks in the darks 
for the UVUDF. In this way more hot pixels are flagged 
than are real, but this ensures that most of the hot pix¬ 
els are masked. The final hot pixel mask is therefore 
equivalent to the green line in Figure im and the aver¬ 
age background flux is shown as the orange line in Figure 

m 

A.5. Impact of Darks On Science Photometry 

In order to test the effects of the blotchiness on the 
catalog photometry, an empty source catalog was gener¬ 
ated in which the source apertures computed from the 
detection image were placed on blank regions in the UV 



Figure 17. Median values of the distribution of measured F275W 
SN measurements and their standard deviation as a function of 
source aperture area. The green squares and dashed line show this 
for science mosaics without a blotchiness correction, and the blue 
circles and solid line show it when including a blotchiness correc¬ 
tion. If the sky noise distribution were uncorrelated and Gaus¬ 
sian, the distribution of the SN would have a median of zero and 
a standard deviation equal to one, independent of aperture area. 
However, an increase in both the median and standard deviation 
of the SN values is observed as the aperture area increases. The 
application of the improved darks significantly improves the pho¬ 
tometry in the NUV, bringing the median closer to zero, and the 
standard deviation closer to one. 
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images. Figure [T7] shows the distribution of the measured 
NUV SN for these blank sources as a function of aperture 
area. If the sky noise distribution were uncorrelated and 
Gaussian, the distribution of the SN for blank sources 
would have a median of zero and a standard deviation 
equal to one, independent of aperture area. However, an 
increase in both the median and standard deviation of 
the SN values is observed as the aperture area increases, 
suggesting that a SExtractor run on the UVUDF images 
would yield measurements with excess flux and underes¬ 
timated uncertainties for large apertures. 

This flux excess and uncertainty underestimation is 
caus ed by the blotchiness pattern described in Section 
IA.3I The number of positive 3cr outliers in the distribu¬ 
tion of flux values in a blank region of one of the UVUDF 
images is ~ 4 times that of the negative 3a outliers. 
Since SExtractor computes the sky background using 
a 3cr clipped mean, these outliers will not be accounted 
for in the sky subtraction and will be falsely attributed 
to flux from the source. In addition, the blotchy pat¬ 
tern introduces a noise pattern that is not accounted for 
in the RMS images used to measured the uncertainty of 
sources. 

This effect is significantly stronger for science mosaics 
using darks without a correction for the blotchiness, com¬ 
pared to mosaics using darks that do. The green squares 
and dashed line in Figure [17] correspond to measurements 
on science data using uncorrected darks, and the blue 
circles and solid line represents the same for darks cor¬ 
rected for the blotchiness. Specifically, the uncorrected 
darks are equivalent to the CDBS darks, except that they 
include a correction for the background gradient and the 
hot pixels. The blotch iness corrected darks are created 
as described in Section lAT4l and include the new average 
dark background per anneal cycle, which removes much 
of the blotchy pattern. After correcting the science mo¬ 
saics for the blotchiness via the darks, the median S/N 
drops from ~ 0.9 to ~ 0.1 for 1000-pixel apertures. Even 
after the correction for blotchiness, there are still some 
residual systematic effects at very large apertures, so 
photometry in that regime should be approached with 
caution. No corrections are applied in the catalogs for 
the remaining blotchy pattern, as only a small number 
of sources have such large areas, and the blotchiness cor¬ 
rected images behave sufficiently well for the majority of 
the sources in the catalog. 
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